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ABSTRACT 


This  dissertation  investigates  power  requirements  for  a  Free  Electron  Laser  to 
bum  through  various  missile  radome  materials.  It  also  includes  computer  simulation 
results  for  several  FEL  system  configurations  designed  to  achieve  maximum  power  while 
maintaining  strict  energy  spread  constraints. 

The  method  used  to  determine  power  requirements  to  bum  through  materials  was 
to  use  the  Thomas  Jefferson  National  Accelerator  Facility’s  Free  Electron  Laser  to 
conduct  material  damage  experiments.  As  the  laser  was  improved  and  increased  in 
power,  the  laser  spot  sizes  on  the  target  materials  were  increased  while  maintaining  a 
constant  irradiance.  The  key  results  from  these  experiments  included  determining 
minimal  spot  sizes  that  can  be  used  for  future  experiments,  and  validation  that  an 
irradiance  level  of  10  kW/cm2  can  bum  through  most  missile  radome  materials  in  a  few 
seconds. 

The  computer  simulations  involved  changing  various  parameters  of  a  FEL  such  as 
electron  energy  levels,  pulse  lengths,  magnetic  field  strengths,  desynchronism,  as  well  as 
several  other  parameters,  to  determine  the  best  possible  configuration  to  achieve  the 
desire  power  levels  and  energy  spread  requirements  for  development  of  a  megawatt  size 
FEL.  The  results  indicate  that  for  the  proposed  designs,  both  the  required  power  and  the 
required  energy  spread  limit  can  be  met. 
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I. 


INTRODUCTION 


Since  the  end  of  the  Cold  War  the  United  States  Navy 
has  adopted  operational  strategies  that  follow  the  ideas 
presented  in  From  the  Sea  and  Forward .  .  .From  the  Sea  [Refs. 
1  and  2]  .  The  central  idea  in  those  documents  is  the 
necessity  of  a  new  maritime  concept,  one  in  which  the  Navy 
plays  a  fundamental  role  in  shore  operations  being  carried 
out  by  our  armed  forces.  As  a  result  Navy  ships, 
submarines,  and  aircraft  can  be  expected  to  operate  in  the 
littoral  regions  of  enemy  territories,  with  the  increased 
potential  for  encountering  hostile  forces.  Even  with  the 
advanced  defensive  systems  on  board  most  of  our  surface 
ships,  the  decreased  reaction  times  available  in  littoral 
operations  makes  them  considerably  more  vulnerable  to  anti¬ 
ship  cruise  missiles  (ASCMs)  and  attacks  by  fast  attack 
boats.  Additionally,  the  close  proximity  of  hostile, 
neutral,  and  friendly  forces  makes  precision  targeting 
extremely  important. 

Building  on  this  new  maritime  concept,  the  Chief  of 
Naval  Operations  recently  released  the  Navy  Strategic 
Planning  Guidance  (April  2000)  which  prioritizes  the 
capabilities  needed  to  support  such  a  maritime  concept.  The 
priority  levels  range  from  Priority  (I) ,  "Those  capabilities 
that  directly  support  or  enhance  the  enduring  core  naval 
competencies  without  which  severe  strategic  risk  would  be 


1 


incurred",  down  to  Priority  (V)  ,  "Those  capabilities  that 
directly  support  or  enhance. . .  without  which  minimal 
strategic  risk  would  be  incurred."  [Ref.  3,  pgs.  40,  61  and 
62] 

Listed  under  priority  (I)  capabilities  are: 

•  The  Fleet-wide,  point  defense  capability  to  achieve 
high-probability  hardkill  against  sub-  and  super¬ 
sonic  cruise  missiles. 

•  The  capability  to  provide  area  air  and  missile 
defense  against  emerging  threats,  including  advanced 
cruise  missiles. 

•  The  capability  to  detect,  identify,  track  and 
destroy  high  numbers  of  small  craft  in  the 
littorals . 


According  to  the  Office  of  Naval  Intelligence  (ONI) , 
the  trends  in  ASCMs  will  include: 

•  Significant  increase  in  missile  terminal  velocity 
from  predominately  sub-sonic  speeds  today  to 
predominately  supersonic  by  2010  and  hypersonic  by 
2020. 

•  Complex  terminal  maneuvers  and  seekers  to  tax  point 
defense  systems. 


Additionally,  ONI  concludes  that  the  capabilities  most 
likely  to  be  encountered  include  massed  small  boat  attack, 
unalerted  single  salvo  cruise  missile  attack,  short-range 
ballistic  missiles,  and  terrorist  attacks.  [Ref.  3,  pgs.  8 
and  14]  Nearly  70  nations  have  deployed  sea  and  land 
launched  ASCMs,  and  20  nations  have  air  launched  ASCMs  [Ref. 
4,  pg.  4] . 
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Super-sonic  and  highly  maneuverable  missiles  are 
already  held  by  known  hostile  nations,  including  the  Russian 
built  SSN-22  Sunburn  which  has  a  cruise  speed/ terminal  speed 
of  mach  2.5/mach  4.5  [Ref.  5,  pg.  121] .  Though  the  terminal 
speed  alone  is  probably  enough  to  greatly  reduce  the 
effectiveness  of  most  self  defense  systems,  this  missile 
also  includes  a  terminal  maneuver.  This  will  reduce  the 
effectiveness  of  self-defense  systems  such  as  the  Navy's 
Close  in  Weapon  System  (CIWS) ,  which  is  especially 
vulnerable  to  these  maneuvers  [Ref.  6,  pg.  179] . 

One  of  the  nations  that  hold  this  missile  is  Iran, 
which  has  the  land  launched  version  [Ref.  5,  pg.  325].  If 
launched  in  the  Straits  of  Hormuz,  which  is  approximately  40 
nautical  miles  wide,  it  would  take  less  than  50  seconds  to 
impact  a  ship  transiting  approximately  20  nautical  miles 
away  through  the  central  section  of  the  straits .  Current 
self-defense  systems  would  have  an  extremely  difficult  time 
engaging  such  a  threat  within  the  50  second  window, 
especially  if  it  was  an  unalerted  launch.  It  is  threats 
such  as  these  that  make  the  idea  of  a  Directed  Energy  Weapon 
(DEW)  stand  out  as  a  possible  solution. 

One  such  DEW  would  be  a  High  Energy  Laser.  In  a  study 
carried  out  by  the  National  Research  Council  to  determine 
Technology  needed  for  the  Navy  and  Marine  Corps  for  the 
years  2000  -  2035,  a  tunable  Free  Electron  Laser  (FEL)  is 
specifically  called  out  as  the  laser  of  choice  for  such  a 
weapon  due  to  its  ability  to  be  tuned  to  wavelengths  that 
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offer  minimal  atmospheric  absorption  and  scattering  [Ref.  6, 
pg.  184] .  A  secondary  advantage  of  using  a  EEL  is  that  its 
gain  medium  is  a  vacuum,  therefore  there  is  no  requirement 
for  the  removal  of  large  amounts  of  waste  heat  from  the 
lasing  medium  as  with  a  solid  state  laser,  or  the  removal  of 
toxic  effluents  associated  with  chemical  lasers. 

Several  topics  are  presented  in  this  dissertation  on 
the  use  of  a  FEL  as  a  Directed  Energy  Weapon.  Chapter  II 
discusses  the  results  of  material  damage  experiments 
conducted  at  Thomas  Jefferson  National  Accelerator  Facility 
(TJNAF)  using  a  FEL.  Some  of  this  work  was  originally 
presented  by  the  author  at  the  21st  International  FEL 
Conference  in  Hamburg,  Germany,  and  is  expected  to  be 
published  in  Ref .  7 . 

Chapter  III  discusses  classical  FEL  Theory.  Starting 
with  the  Lorentz  force  equation  and  Maxwell's  wave 
equation,  the  dimensionless  FEL  pendulum  and  wave  equations 
are  developed  [Ref.  29].  The  work  in  this  thesis  includes 
modifications  of  the  FEL  interactions  and  gain  to  include 
short  pulse  effects  and  tapering  are  also  developed, 
including  a  modified  pendulum  equation  to  include  the 
effects  of  tapering  of  the  undulator.  This  work  was 
presented  by  the  author  at  the  21st  International  FEL 
Conference  in  Hamburg,  Germany  and  was  published  in  Ref.  8. 
More  recent  results  from  simulations  conducted  by  the  author 
were  presented  at  the  22nd  International  FEL  Conference  in 
Raleigh,  NC,  and  will  be  published  in  Ref.  7. 
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Chapter  IV  includes  the  results  of  computer  simulations 
for  proposed  output  power  upgrades  to  the  TJNAF  FEL.  These 
upgrades  are  the  next  step  in  a  proposed  upgrade  to  1  MW 
output  power.  This  work  was  analyzed  by  the  author  and  was 
presented  at  the  20th  International  FEL  Conference  in 
Williamsburg,  VA,  and  22nd  International  FEL  Conference  in 
Raleigh,  NC.  Results  from  20  kW  simulations  were  published 
Ref.  .9,  and  results  for  the  10  kW  simulations  will  be 
published  in  Ref.  7. 

The  last  chapter  discusses  the  effect  of  tapering  the 
undulator ,  and  specifically  calls  out  the  differences  and 
advantages  of  positive  vs.  negative  taper.  Simulations 
conducted  by  the  author  resulted  in  surprising  positive 
results,  and  were  presented  by  him  at  the  21st  and  22nd 
International  Free  Electron  Laser  Conferences,  and  published 
in  Refs.  7  and  8. 
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II.  MATERIAL  DAMAGE  EXPERIMENTS 


A.  THEORY 

The  interaction  between  lasers  and  matter  is  a 
complicated  issue  generally  described  by  non-linear  and 
unpredictable  effects.  The  basic  kill  mechanism  for 
destroying  a  missile  is  to  cause  a  catastrophic  failure  of 
one  of  the  central  components  such  as  the  warhead,  control 
and  navigation  electronics,  or  the  high  pressure  chamber  in 
the  propulsion  section,  or  to  create  aerodynamic  instability 
beyond  the  design  limits  of  the  missiles  flight  control 
surfaces.  [Ref.  10] 

The  laser  beam  has  the  unique  ability  to  deliver  very 
high  power  per  unit  area.  As  the  optical  energy  strikes  the 
surface,  it  is  partially  absorbed  and  partially  reflected. 
For  metal  surfaces  and  metal  alloys,  the  reflectivity  is 
typically  very  high  (92%  or  greater)  at  room  temperature. 
However,  for  metals  interacting  with  infrared  wavelength 
radiation,  there  is  a  decrease  in  reflectivity  with 
temperature,  with  an  abrupt  drop  in  reflectivity  as  the 
material  changes  from  solid  to  liquid  state  as  shown  in 
Table  1.  TSL  is  the  melting  temperature  of  the  metal.  [Ref. 
11.]  This  loss  of  reflectivity  causes  a  large  increase  in 
the  amount  of  energy  absorbed,  which  results  in  an  even 
faster  increase  in  temperature.  Provided  the  amount  of 
energy  added  is  greater  than  the  rate  of  energy  diffusion  to 
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surrounding  cooler  areas,  the  surface  will  become  hotter  and 
hotter  and  change  from  solid  to  liquid  and  finally  become 
vaporized. 


Table  1.  Absorption  Percentage  vs.  Temp  for  Metals.  (After  Ref.  11) 


Metal 

tb 

(K) 

Absorption 

at  300  K 

Absorption 

at  T 

m 

(Solid) 

Absorption 

at  T 

m 

(Liquid) 

Absorption 

at 

T  +  500K 

xn 

Ag 

1234 

0.4% 

2.0% 

4.3% 

5.4% 

Al 

933 

1.1% 

4.4% 

10.0% 

13.0% 

Au 

1336 

0.6% 

3.5% 

8.3% 

10.3% 

Cu 

1356 

0.5% 

2.8% 

5.8% 

7.3% 

For  typical  radome  materials  there  are  basically  two 
mechanisms  of  destruction.  The  first  is  penetration  by 
melting  and  ablating  the  material.  The  second  is  to  create 
fractures  by  thermally  stressing  the  material. 
Additionally,  if  the  missile  housed  an  infrared  seeker,  the 
radome  would  naturally  be  transparent  in  the  infrared.  In 
that  case,  destruction  of  the  radome  would  not  be  necessary. 
Instead,  the  infrared  laser  would  penetrate  the  radome  and 
interact  directly  with  the  guidance  system  and  other  parts 
behind  the  radome.  [Ref.  12] 
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Some  of  the  physical  processes  that  govern  the  damage 
caused  to  the  material  by  the  laser  include  absorption  and 
reflection,  density  of  the  material,  heat  capacity,  and 
thermal  conduction.  Furthermore,  the  characteristics  of  the 
laser  such  as  irradiation  wavelength,  power  density,  peak 
power,  average  power,  and  pulse  characteristics  play  an 
important  role  as  well. 

1.  Required  Power  to  Destroy  Target 

One  of  the  major  reasons  behind  the  experiments 
discussed  in  this  chapter  was  to  determine  just  how  much 
power  from  a  short-pulse  FEL  is  needed  to  destroy  a  missile 
in  the  few  seconds  allowed  for  an  engagement. 

One  can  estimate  the  required  irradiance  by  assuming 
that  the  laser  burns  through  the  material  by  breaking  the 
cohesive  bonds  of  individual  atoms  allowing  the  removal  of 
the  atoms.  In  actuality,  the  matter  would  most  likely 
disintegrate  in  "chunks"  of  atoms.  When  the  laser  beam 
strikes  the  missile,  not  all  the  energy  is  absorbed,  and 
even  some  of  the  absorbed  energy  makes  its  way  to  cooler 
parts  of  the  material  prior  to  melting.  Additionally,  some 
of  the  power  may  be  wasted  by  heating  some  of  the  atoms 
beyond  the  temperature  needed  to  remove  them. 

One  of  the  materials  irradiated  during  the  experiments 
was  aluminum.  Metals  such  as  aluminum  could  be  expected  to 
make  up  the  body  casing  of  a  cruise  missile.  The  binding 
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energy  of  aluminum  is  approximately  3.4  eV/atom  [Ref.  13, 
pg.  74]  .  If  the  casing  is  assumed  to  be  3  layers  thick, 
with  each  layer  approximately  1  cm  thick,  and  the  laser  spot 
size  on  the  target  is  100  cm2,  the  volume  of  material  to  be 
removed  is  approximately  300  cm3.  The  atomic  spacing  for 
aluminum  is  approximately  2.5  angstroms  [Ref.  13.  pg.  98], 
or  2.5  x  10'8  cm,  which  corresponds  to  the  total  number  of 
atoms  to  be  removed 


_ 300  cm3 _ 

(  2.5  •  10~8cm)3  /  atom 


2  •  1025  atoms. 


(1.) 


The  total  energy  required  is  therefore  estimated  at 


E  =  (2  •  1025 atoms)  (3 .4eV  /  atom)  (l .  6  •  10",9J  /  ev) 

\Z  *  ) 

~  10  MJ 

To  deliver  10  MJ  of  energy  in  3  to  4  seconds  requires 
approximately  3  MW  at  the  target,  assuming  all  the  energy  is 
absorbed.  Three  MW  divided  over  100  cm2  gives  a  required 
irradiance  of  0  =  30  kW/cm2.  The  actual  required  irradiance 

level  would  be  determined  by  the  amount  of  reflected  energy. 

A  second  method  to  estimate  the  required  irradiance  is 
to  determine  the  amount  of  energy  needed  to  bring  the 
material  to  its  vaporization  temperature.  If  that  the 
energy  is  delivered  at  a  rate  much  greater  than  the  heat 
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loss  through  diffusion,  the  required  energy  can  be 
determined  from 

<Do  =  pz  (c[Tm  -  T0]  +  AHra  +  C[TV  -  Tm]  +  AHv)  (3  . ) 

where  <I>o  is  the  required  flux  density,  p  =  mass  density,  z  = 
thickness  of  material  to  be  burned  through,  C  =  Specific 
Heat,  T  =  melting  temperature,  T  =  ambient  temperature,  Tv 
=  vaporization  temperature,  AHm  =  latent  heat  of  melting, 
and  AHv  =  latent  heat  of  vaporization  [REF.  14,  pg .  167]. 

For  aluminum  the  specific  values  are  p=2700  kg/m  ,  z  -  3 
cm,  C  =  896  J/kg-K,  Tm  =  855  K,  T0  =  300  K,  Tv  =2750  K,  AHm  = 
4  x  10b  J/kg,  and  AHv  =  10.8  x  106  J/kg.  For  this  rough 

estimate  the  same  value  for  specific  heat  has  been  used  for 
both  the  solid  and  liquid  phases,  although  in  actuallity 
these  numbers  may  be  different.  Using  these  values  in 
equation  (3.)  gives  a  required  flux  density  of 

$  =  1  GJ  /  m2  (4  • ) 

O 

or 

<Do  =  100  kJ  /  cm2  (5.) 

This  much  energy  delivered  over  a  3  second  engagement 
requires  an  irradiance  of  approximately  <J>0=35  kW/cm  , 
corresponding  to  a  total  beam  power  of  3 . 5  MW  in  a  100  cm2 
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spot  on  the  target.  This  result  is  consistent  with  the  3  MW 
estimate  developed  with  the  first  approximation  method. 

A  second  material  irradiated  during  the  experiments  was 
Slip-Cast  Fused  Silica,  a  furnace  tile  material.  The 
specific  parameters  are  p-  2200  kg/m3,  z  =  3  cm,  C  =  920 
J/kg-K,  Tn  =  1980  K,  T0  =  300  K,  Tv  =  2200  K,  AHm  =  1.5  x  106 
J/kg,  AHv  =  2.2  x  106  J/kg.  Using  these  numbers  in  equation 

(3.)  shows  that  35  kJ/cm2,  corresponding  to  10  kW/cm2  over  a 
3  to  4  second  irradiation  in  order  to  burn  through  the 
material.  This  is  three  times  smaller  than  for  aluminum 
because  the  latent  heat  of  vaporization  for  aluminum  is  five 
times  that  of  the  silica,  and  the  vaporization  temperature 
of  aluminum  is  also  higher.  The  thickness  of  the  Slip-Cast 
Fused  Silica  used  in  the  experiments  was  actually  only  9  mm, 
about  1/3  of  the  thickness  of  the  aluminum  used  in  this 
example.  Calculating  the  required  energy  for  the  9  mm 
sample  shows  that  the  total  energy  needed  is  only  11  kJ/cm2, 
corresponding  to  a  required  irradiance  of  approximately 
<J>0  =  3  kW/cm2  for  a  3  to  4  second  engagement. 

These  approximations  should  provide  a  good  estimate  of 
the  energy  that  is  actually  needed  to  melt  through  the 
missile.  As  reported  later  in  this  chapter,  an  irradiance 
level  of  10  kW/cm2  was  enough  to  melt  through  an  aluminum 
sample  that  was  cut  to  the  correct  size  in  order  to  control 
the  thermal  diffusion.  However,  when  a  sample  of  Slip-Cast 
Fused  Silica  was  irradiated  at  the  same  irradiance  level. 
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burn  through  was  only  achieved  after  an  extended  time  of 
almost  two  minutes.  A  possible  explanation  for  this  will  be 
discussed  later  in  this  chapter. 

Experiments  conducted  in  the  1970 's  and  1980 's  in 
conjunction  with  the  MIRACL  program  indicate  that  a  power 
density  in  the  tens  of  kilowatts  is  needed  to  destroy  a 
missile  with  a  dwell  time  of  a  few  seconds.  Experiments 
conducted  by  the  Air  Force  during  the  same  time  period 
determined  that  a  flux  of  50  kW/cm2  was  needed  to  damage  the 
missile  structure  in  less  than  one  second  [Ref.  15]  .  Both 
of  these  sets  of  tests  seem  to  validate  the  order  of 
magnitude  approximations  made  in  this  section.  For  this 
reason  most  of  the  experiments  discussed  later  in  the  report 
were  conducted  at  an  irradiance  level  of  approximately 
10  kW/cm2. 


2.  Scaling  Laws  and  Minimal  Spot  Size 

Since  there  is  no  MW-class  FEL  to  perform  full-scale 
experiments,  scaling  is  the  only  way  to  determine  the 
effectiveness  of  a  FEL  weapon.  As  discussed  in  the  previous 
section,  a  power  density  of  approximately  <1>  =  10  kW/cm2  over 
a  spot  size  A  =  100  cm2  is  required  to  destroy  a  missile 
with  a  dwell  time  of  a  few  seconds.  Scaling  laws  would 
allow  predictions  of  large  area  damage  from  small  area 
experiments  with  the  confidence  that  the  two  were  correlated 
correctly. 
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To  achieve  a  power  density  of  O  =  10  kW/cm2,  a  100  W 
FEL  must  use  a  spot  size  of  1  mm2,  while  a  1  kW  laser  uses  a 
larger  spot  size  of  10  mm2.  Scaling  of  the  laser  damage 
will  only  work  if  the  thermal  diffusion  is  independent  of 
spot  size,  which  is  only  true  for  laser  spots  larger  than 
the  characteristic  length  for  thermal  diffusion.  Schriempf 
[Ref.  16]  calculates  the  characteristic  thermal  diffusion 
length  to  be 


Dl  =  2jxt~e 


(6.) 


where 


is  the  thermal  diffusivity  and 


7t_  K2AT2 


(7.) 


where  te  is  the  time  required  to  bring  the  material  from 
ambient  temperature  to  melting  temperature,  K  is  the  thermal 
conductivity,  AT  is  the  temperature  change  required  to 
reach  melting  temperature,  and  <I>0  is  the  energy  flux  being 
absorbed  by  the  material.  The  thermal  diffusion  length  DL 
represents  the  distance  required  for  temperature  T  to  drop 
to  1/e  times  its  central  value.  In  the  semi-infinite 
approximation  used  by  Schriempf,  radial  heat  flow  is 
ignored.  In  order  for  this  to  be  valid,  the  spot  size  must 
be  much  larger  than  DL,  or  the  target  diameter  smaller  than 
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If  these  conditions  are  not  met,  heat 


Dl,  diameter  <  DL. 
will  diffuse  away  from  the  laser  spot,  and  the  area  will  not 
be  heated  effectively. 


3 .  Pulse  Train 

FELs  produce  short,  powerful  pulses  with  a  rapid 
repetition  rate  in  the  MHz  range.  This  can  be  compared  with 
the  typical  high  peak  power  short  pulse  lasers  that  operate 
at  very  low  repetition  rates  of  Hz  or  kHz . 

When  the  laser  power  is  delivered  in  pulses,  the  peak 
power  increases  and  may  cause  new  effects  beyond  thermal 
heating.  Short  pulsed  lasers  (microseconds  and  shorter) , 
and  ultra-short  pulsed  lasers  (picoseconds  and  shorter)  are 
capable  or  producing  impulse  damage  in  addition  to  the 
thermal  damage.  The  high  peak  power  can  cause  rapid 
vaporization  at  the  target  surface,  and  subsequent  recoil 
from  the  vapor  blowoff  can  form  a  strong  pressure  wave.  The 
pressure  wave  forms  a  reflected  wave  at  the  backside  of  the 
material,  and  the  superposition  of  the  two  waves  can 
overstress  the  material  and  cause  catastrophic  failing  of 
the  material.  [Ref.  17] 

Another  effect  occurs  where  an  ultra-short  laser  pulse 
delivers  energy  to  a  metal  at  such  a  fast  rate  that  the 
metal  lattice  cannot  respond,  but  the  electrons  can.  The 
electrons  rapidly  increase  in  temperature  so  that  the 
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difference  between  electron  and  lattice  temperatures  can  be 
as  much  as  a  few  thousand  degrees.  [Ref.  18] 

Experiments  have  shown  that  the  damage  threshold 
fluence,  Eth,  of  a  material  gets  smaller  as  the  pulse  length 
gets  smaller.  Experiments  show  that  there  may  be  as  much  as 
a  factor  of  10  advantage  when  using  shorter  picosecond 
pulses  over  the  larger  100  nanosecond  pulses.  [Ref.  19] 

Research  has  also  examined  the  threshold  for  ultra- 
short  pulses  damaging  dielectrics.  The  same  trend  as  was 
observed  with  metals  was  observed  in  the  dielectrics.  For 
the  longer  pulses,  the  damage  threshold  decreases  with 
decreasing  pulse  length  tp  as  tp1/2,  down  to  very  short, 
picosecond  pulses.  [Refs.  20,  21] 

In  the  experiments  conducted  at  TJTSJAF,  the  higher 
fluence  pulses  contained  0.12  J/cm2  of  energy,  which  is 
slightly  higher  than  the  ablation  threshold  for  picosecond 
pulses  on  metals  [Ref.  22]  .  Some  ablation  of  the  material 
may  have  taken  place.  However,  the  fluence  level  was  not 
high  enough  to  create  the  impulse  damage  and  pressure  waves 
discussed  in  the  paragraphs  above. 

B.  MATERIALS  USED  IN  EXPERIMENTS 

In  addition  to  the  aluminum  samples  tested  in  the 
experiments,  five  other  materials  were  also  used.  The 
materials  chosen  were  either  actual  missile  nose  cone 
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materials,  or  materials  that  were  very  similar  to  missile 
nose  cones.  The  materials  are  [Refs.  23,  24]: 

1.  Phenolic  Resin: 

Radome  material  from  a  Standard  ARM  Missile. 

2.  Pyroceram: 

A  furnace  coating  type  material  used  in  SM-1 
and  SM-2  missile  radomes. 

3.  Slip-Cast  Fused  Silica: 

Commercially  available  furnace  material. 
Similar  to  material  used  as  a  radome  for 
Patriot  missile  system. 

4 .  Polyimide  Fiberglass : 

A  high  temperature  fiberglass  used  on 

supersonic  missile  radomes. 

5 .  F2  Epoxy : 

Five  plies  of  fiberglass  cloth  with  an  epoxy 
binder.  This  type  of  material  is  used  in  the 
Soviet  STYX  type  missiles. 

1.  Thermal  Diffusion  Calculations 

The  first  step  in  preparing  for  the  experiments  was  to 
calculate  the  thermal  diffusion  length  for  each  material. 
The  first  calculation  and  test  was  made  for  a  small  piece  of 
aluminum  (Al-6061)  .  This  material  had  density  p  =  2700 
kg/m3,  specific  heat  C  =  896  J/kg-K,  thermal  conductivity  K 
=  180  W/m-K,  thermal  diffusivity  K  =  7.44  x  10'5  m2/s, 


17 


melting  temperature  Tm=  855K,  and  power  density  <3>0  =  10s  W/m2 
=  10  kW/cm2.  From  equation  (6.),  the  thermal  diffusion 

length  is  D  =  2  mm,  with  a  time  t=0.01  seconds  required  to 
begin  the  melting  process.  In  order  to  melt  through  an 
aluminum  sample,  the  laser  spot  must  have  an  area  greater 
than  n( 2  mm)2  =  10  mm2.  This  would  make  the  spot  large 
enough  so  that  the  heat  arriving  at  the  center  of  the  spot 
would  not  be  able  to  diffuse  to  the  cooler  regions  outside 
the  spot  before  melting  could  occur.  Alternately,  the 
target  itself  can  be  made  small  with  the  diameter  of  the 
sample  diameter  ~  D  so  as  to  prevent  the  heat  from  escaping 
to  cooler  regions  of  material  before  heating  could  take 
place . 

These  calculations  were  experimentally  verified  with 
samples  of  Al-6061.  A  sample  with  a  1  cm  diameter  was 
irradiated  with  a  laser  spot  size  of  1  mm2  (much  smaller 
than  the  required  10  mm2)  without  melting  after  several 
minutes.  Other  samples  with  a  diameter  of  2  mm  (very  close 
to  the  laser  spot  size  itself)  were  melted  in  a  few  seconds. 

Another  target  irradiated  was  Slip-cast  Fused  Silica 
(Sio2)  .  A  calculation  of  the  thermal  diffusion  length 
associated  with  heating  the  sample  to  its  melting 
temperature  of  1980  K  was  performed  using  values  p  —  2200 
kg/m3,  C  =  920  J/kg-K,  K  =  1.26  W/m-K,  K  =  5  x  10'7  m2/s,  Tn= 
1980K,  <I>0  =  10s  W/m2  =  10  kW/cm2.  Since  this  material  has  a 
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low  thermal  dif fusivity,  its  thermal  diffusion  length  is 
also  small,  DL=0.02  mm.  An  insulating  material  such  as 
fused  silica,  the  scaled  laser  spot  can  have  an  area  greater 
than  71(0.02  mm)2  =  0.001  mm2,  which  was  easily  attained  with 
the  1  mm2  or  larger  beams  used  at  TJNAF. 

2.  Experimental  Results 

The  experiments  conducted  from  March  1999  through  March 
2000  are  reported  in  great  detail  by  McGinnis,  et  al.  [Ref. 
25]  Results  from  two  of  the  materials  are  reviewed  in  this 
section. 


a)  Slip-Cast  Fused  Silica 

One  of  the  more  thorough  tests  was  conducted  on 
Slip-Cast  Fused  Silica,  Si02.  The  front  of  the  sample  after 
all  irradiations  is  shown  in  Figure  1,  with  the  numerical 
results  recorded  in  Table  2 . 
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Figure  1.  Slip-Cast  Fused  Silica 


Table  2.  Irradiation  Results  of  Slip-Cast  Fused  Silica. 


Run 

Number 

Wavelength 

(pm) 

PRF 

(MHz) 

Airflow 

(mph) 

Average 

Irradiance 

(kW/cm2) 

Laser 

Beam 

Radius 

(mm) 

Penetration 

Rate 

(mm/s) 

1 

3.1 

37.4 

83-86 

10 

1.2 

1.32 

2 

3.1 

37.4 

NO 

10 

1.2 

7.5 

3 

4.8 

37.4 

NO 

9 

0.87 

0.2 

4 

4.8 

37.4 

NO 

9 

0.87 

0.081 

5 

4.8 

37.4 

NO 

9 

0.87 

0.2 

6 

4.8 

37.4 

NO 

9 

0.87 

0.125 

7 

4.8 

37.4 

NO 

9 

0.87 

0.081 

8 

4.8 

37.4 

NO 

9 

0.87 

0.35 

9 

4.8 

37.4 

NO 

500 

0.08 

3.0 

10 

4.8 

37.4 

NO 

500 

0.08 

3.0 
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As  can  be  seen  in  Table  2,  there  is  quite  a 
variation  in  the  penetration  rates  for  the  different  runs. 
As  the  exposure  time  was  increased,  more  smoke  and  melted 
material  filled  the  hole  and  blocked  the  path  of  the  laser 
beam  causing  the  penetration  rate  to  decline  over  time. 
Apparently  9  kW/cm2  is  not  a  high  enough  irradiance  for  this 
material  to  reach  vaporization  temperature.  The  heat  was 
able  to  diffuse  away  quickly  enough  that  the  material 
reached  a  steady  state  in  the  liquid  form,  and  would  not 
clear  for  additional  penetration.  However,  runs  9  and  10 
burned  through  very  quickly  indicating  that  500  kW/cm2  is 
more  than  adequate  to  bring  the  material  to  vaporization 
temperature . 

Since  run  four  was  the  only  run  to  burn  completely 
through  the  material  at  a  power  density  of  9  kW/cm2,  it  is 
worth  a  closer  look.  A  digital  picture  of  damage  from  run 
four  was  taken  through  an  optical  microscope  as  shown  in 
Figure  2 . 
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Figure  2.  Damage  to  Slip-Cast  Fused  Silica  Run  4. 

The  laser  beam  diameter  narrowed  from  the  front 
surface  to  the  back  surface  of  the  hole;  its  dimensions  were 
1.1  mm  at  the  front  and  approximately  0.7  mm  at  the  back. 
The  actual  hole  (the  section  of  the  sample  with  material  had 
been  completely  removed  during  lasing) ,  was  very  close  to 
those  dimensions  ranging  from  1.06  mm  at  the  front  down  to 
0.7  mm  at  the  back.  However,  damage  to  the  material  was 
much  more  wider,  with  the  melted  portion  at  the  surface  of 
the  sample  measuring  5.0  mm  in  diameter,  and  the  hole  at  the 
back  edge  measuring  2.0  mm  in  diameter. 

Examination  of  the  hole  from  run  four  through  an 
optical  microscope  revealed  a  1  mm  thick  layer  of  melted  and 
rehardened  Si02  filling  the  hole  at  the  back  of  the  sample. 
It  was  clear  from  the  video  and  from  the  rear  power  meter 
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that  burn- through  occurred  in  run  four,  but  melted  material 
solidified  and  resealed  the  hole  at  the  back  of  the  sample. 

The  volume  of  the  hole  in  run  four  is  estimated 

by, 


/•9mm  0 

V  =  Jo  ;zR2(z)dz,  (8.) 

where  the  radius  closely  matches  the  laser  beam  radius  and 
approximately  changes  linearly  as 

R ( z )  =  0.53  mm  -  0.02  z  (9.) 

where  z  is  the  distance  from  the  front  of  the  sample.  Which 
gives  a  volume  of  V  =  5.6  mm3  or  5.6  x  10’3  cm3  .  By  a 
similar  calculation,  the  volume  of  the  entire  damaged 
region,  including  the  melted  and  rehardened  portion,  is 
estimated  to  be  V  =  92  mm3.  Based  on  the  density  of  fused 
silica,  p  -  2.2  grams/cm2,  the  amount  of  material  removed 
was  0.012  grams,  while  the  amount  of  material  damaged 
(melted)  was  0.20  grams.  The  heat  energy  deposited  during 
run  four  is  given  by 


E  =  Pt  =  OAt,  (10. ) 

where  A  =  71(0.53  mm)2,  0  =  9  kW/cm2,  and  the  irradiation 

time  t  =  110  s,  so  that  E  =  9.7  kJ  deposited. 

The  predictions  made  in  section  A.l.  of  this 
chapter  indicated  that  this  material  could  be  burned  through 
in  about  3  seconds  at  3  kW/cm2,  or  just  over  1  second  at  the 
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9  kW/cm2  used  in  these  experiments.  Using  those  numbers  in 
equation  (10.)  gives  an  estimation  of  the  energy  that  would 
need  to  be  deposited  to  melt  through  this  sample  is 

E  =  Irradiance  •  Beam  Area  •  time 

=  (3  kW  /  cm2)  •  7t  (5.5  •  10"2cm)2  -(3  s)  (11. ) 

=  0.95  kJ 


vice  the  9.7  kJ  that  was  actually  deposited. 

The  amount  of  energy  needed  can  also  be  calculated 
using  the  first  method  in  section  A.l.  That  is  by 
calculating  the  energy  needed  to  break  each  atomic  bond  in 
the  volume  of  the  Si02  affected.  The  atomic  weight  of  Si02 
is  60  amu,  or  1  x  10"22  grams.  Since  the  density  of  Si02  is 
2.2  grams/cm3,  the  number  of  molecules  per  cubic  centimeter 
is  approximately  2  x  1022.  The  required  energy  to  break  all 
the  bonds  assuming  approximately  4  eV/bond  is  estimated  as 

E  =  (  2.2  •  1022  bonds  /  cm3)  (  4eV  /  bond)  ( 1 .6  •  10-l9J  /  ev) 

=  12.8  kJ  /  cm3. 

From  the  volume  calculated  from  equation  (8.),  the 
total  energy  required  is  predicted  to  be 

E  =  ( 12.8  kJ  /  cm3)(  5.6  •  10~3cm3)  ^ 

«  0.72  kJ 
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which  is  very  close  to  the  0.95  kJ  predicted  in  equation 
(11.),  but  far  from  the  9.7  kJ  delivered  during  the 
irradiation  of  run  4 . 

The  difference  in  the  predicted  and  actual  energy 
requirements  may  have  come  about  as  a  result  of  the  change 
in  state  of  the  material.  The  calculations  used  in  section 
A. 1 .  of  this  chapter  and  in  the  paragraphs  above  assume  that 
the  thermal  diffusion  length  is  small  compared  to  the  beam 
diameter.  This  assumption  is  true  for  Slip-Cast  Fused 
Silica  in  its  normal  state,  where  the  thermal  diffusion 
length  was  calculated  in  section  B.l.  to  be  only  0.02  mm. 
However,  it  is  possible  that  once  the  material  melted,  the 
thermal  conductivity  increased  dramatically,  increasing  the 
thermal  diffusion  length  substantially.  This  increase  in 
thermal  conductivity  would  result  in  a  dramatic  increase  in 
the  rate  of  heat  diffused  away  from  the  laser  spot,  and 
prevent  the  material  from  reaching  the  vaporization 
temperature.  This  increase  seems  to  be  what  is  taking  place 
since  even  after  being  irradiated  for  110  seconds,  the 
material  still  contains  a  large  volume  of  melted  and  re¬ 
solidified  material. 

This  postulate  can  be  checked  by  calculating  the 
energy  deposited  during  run  9,  which  was  at  a  much  higher 
irradiance  and  succeeded  in  melting  through  the  hole  in  less 
than  3  seconds.  The  amount  of  energy  deposited  for  this  run 
was 
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(14.) 


E  =  Ir radiance  •  Beam  Area • time 

=  (  500  kW  /  cm2)  •  #(8.0  •  10_3cm)  -(3s) 

=  0 . 3  k  J, 

again  showing  an  amount  of  energy  much  smaller  than  the  = 
10  kJ  deposited  during  burn  through  for  run  4.  This 
reinforces  the  postulate  that  the  energy  must  arrive  at  a 
rate  high  enough  to  raise  the  temperature  of  the  material  to 
vaporization  temperature  before  the  molten  material  can  re- 
solidify.  It  appears  that  almost  the  entire  9.7  kJ  of 
energy  delivered  during  run  4  was  wasted  because  it  was 
delivered  too  slowly. 

In  a  realistic  directed  energy  application,  the 
irradiance  used  in  run  9  would  not  be  needed.  A  much  larger 
spot  size  of  100  cm2  would  not  allow  for  significant  thermal 
diffusion  away  from  the  damaged  area. 

The  final  experiments  were  conducted  in  March 
2000.  The  purpose  of  the  March  2000  experiments  was  to 
increase  the  size  of  the  laser  beam  spot  while  maintaining 
the  same  irradiance  of  approximately  10  kW/cm2.  To  do  this 
the  average  beam  power  was  increased  to  500  Watts.  Air  flow 
was  also  added  to  these  experiments  to  determine  if  airflow 
across  the  material  would  effect  the  damage  process. 

Burn-through  of  the  material  occurred  for  both  of 
the  runs  conducted  to  the  Slip-Cast  Fused  Silica  during  this 
experiment.  The  burn-through  time  for  run  1  was  6.8  sec  and 
for  run  2  was  1.2  sec.  In  run  1,  the  damage  diameter  on  the 
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front  face  of  the  material  sample  was  4 . 4  mm,  and  on  the 
back  face  was  3.1  mm.  In  run  2  the  damage  diameter  on  the 
front  face  was  5.6  mm  and  on  the  back  face  1.8  mm.  The 
effect  of  airflow  was  to  reduce  the  front  face  damage 
diameter  and  to  increase  the  burn- through  time. 

The  volume  of  the  entire  damaged  region  in  run  1 
is  estimated  by. 


V,T  =  71  •  RfT(z)  •  dz  ,  (15  . ) 

where  z  is  the  depth  of  the  hole.  The  radius  changes 
approximately  linearly  as 

R1t(z)  =  2.2mm  -  0.07  •  z  (16.) 

The  volume  of  the  entire  damaged  region  in  run  2 
is  also  estimated  by  equation  (15.),  where  the  radius 
changes  approximately  linearly  as 

R2t(z)  =  2.8mm  -  0.21  •  z  (17.) 

After  doing  the  above  calculations,  the  volume  of 
the  entire  damage  region  in  run  1  is  V1T=100  mm3  and  in  run 
2  is  V2t=105  mm3.  The  density  of  the  fused  silica  is  p=  2.2 
grams /cm3,  therefore  mass  of  the  entire  damaged  region  in 
run  1  is  mlT=  0.22  grams  and  in  run  2  is  m2T=0 .23  grams. 


27 


The  volume  of  the  holes  in  runs  1  and  2  are 
estimated  using  radii 


R1h(z)  =  1.6mm  -  0.17  •  z 


(18.) 


and 


R^z)  =  0.75mm  -  0.04  •  z  (19.) 

The  calculated  volumes  of  the  hole  in  run  1  is 
V1h=35  mm3  and  in  run  2  is  V2H=8.9  mm3.  The  mass  of  the 
material  removed  creating  a  hole  in  run  1  is  m1H=0.08  grams 
and  in  run  2  is  m2H=0.02  grams. 

These  calculations  show  that  the  addition  of  air 
flow  had  a  negligible  effect  on  the  overall  damaged 
regions,  since  both  have  approximately  the  same  volume. 
Air  flow  was  substantial  in  affecting  the  amount  of  material 
removed  since  the  hole  volume  increased  from  V2H=8.9  mm3  to 
V1h=35  mm3  when  the  air  flow  was  blowing  across  the  material. 
However,  the  air  flow  also  decreased  the  burn  through  rate. 
So,  though  the  hole  was  bigger,  it  took  longer  to  create. 

Comparing  the  results  of  the  runs  1  and  2  from  the 
March  2000  experiments  with  those  of  runs  3  through  8  in 
Table  2,  it  can  be  seen  that  shifting  the  wavelength  from 
A  =  4.8  Mm  to  A  =  3.1  M™'  and  increasing  the  average  power 

from  100  Watts  to  500  Watts  had  a  significant  effect  on  the 
penetration  rates.  Those  penetration  rates  increased 
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dramatically  from  0.081  mm/s  to  1.32  mm/s  with  the  presence 
of  airflow  and  to  7.5  mm/s  without  it. 

A  possible  explanation  for  the  success  of  the 
March  2000  runs  is  the  increase  in  beam  radius  from  0.087cm 
to  0.12cm.  This  increase  in  beam  area  could  be  enough  so 
that  even  with  the  larger  thermal  diffusion  length  estimated 
in  the  preceding  paragraphs,  the  new  beam  area  was  large 
enough  to  ensure  enough  heat  was  retained  to  reach 
vaporization. 


b)  Polyimide  Fiberglass 

Experiments  were  also  conducted  using  Polyimide 
Fiberglass.  The  front  and  back  of  the  sample  after  all 
experiments  were  conducted  are  shown  in  Figure  3  and  Figure 
4.  Table  3  shows  the  data  for  all  the  experiments  conducted 
on  this  material  over  the  year  long  period.  Note  that  runs 
10  through  15  were  conducted  at  500  W/cm2,  while  all  the 
other  runs  were  conducted  at  10  kW/cm2. 

In  Figure  4,  we  see  the  backside  of  the  sample 

showing  that  all  of  the  10  kW  /  cm2  irradiations  completely 
penetrated  the  sample.  (Note:  Hole  3  in  Figure  3  has 
mistakenly  been  numbered  2  in  Figure  4)  .  As  can  also  be 

seen,  none  of  the  500  W/cm2  irradiance  runs  (irradiations 
10  to  15  of  Figure  3)  penetrated  the  sample. 
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Figure  3.  Polyimide  Fiberglass  (Front  View). 
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Figure  4.  Polyimide  Fiberglass  (Back  View). 
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Table  3.  Polyimide  Fiberglass  Irradiation  Results. 
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During  runs  1  through  3 ,  only  run  3  actually 
burned  through  the  sample.  These  three  runs  have  the  same 
irradiation  parameters  with  those  of  runs  7,8,  and  9  of 
Table  3  except  for  the  PRF  and  the  wavelength.  Comparing 
these  runs,  it  appears  that  the  shorter  (A,  =  3.1pm) 
wavelength  combined  with  the  lower  PRF  is  much  more 
effective,  resulting  in  500%  higher  penetration  depth  rate. 
The  lower  PRF  may  be  more  efficient  because  it  has  more 
energy  per  pulse  for  the  same  pulse  length. 

For  the  holes  that  burned  through,  the  presence  of 
airflow  resulted  in  larger  holes  and  an  increased 
penetration  depth  rate  than  those  attained  without  air . 

The  larger  spot  size  irradiations  conducted  with 
500  W  average  power  (runs  16  to  21),  caused  3.5  times  higher 
penetration  depth  rate  than  the  ones  using  100  W  (runs  4  to 
9),  even  though  both  were  at  the  same  irradiance  level,  PRF, 
and  wavelength.  This  could  be  because  the  larger  spot 
allows  material  to  exit  the  region  quicker,  or  because  the 
larger  spot  better  meets  the  thermal  diffusion  requirement. 
Investigation  with  a  microscope  reveals  that  there  is  no 
evidence  of  melted  or  rehardened  material  inside  the  holes. 


C.  CONCLUSIONS 

This  chapter  describes  some  of  the  first  measurements 
of  laser  damage  on  missile  materials  from  a  newly  developed 
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TJNAF  Free  Electron  Laser.  A  much  more  detailed  report  of 
all  the  experiments  conducted  over  a  year  long  period  has 
also  been  prepared  [Ref.  25]. 

During  the  experiments  conducted  on  March  12th  and  March 
23rd,  1999,  the  only  parameter  changed  was  the  pulse 
repetition  frequency,  which  caused  the  pulses  used  in  the 
later  experiments  to  have  twice  the  fluence  per  pulse  as  the 
pulses  on  the  12th.  The  burn  through  rate  for  the 
experiments  of  March  23rd  was  faster,  and  it  is  postulated 
that  the  increase  occurred  because  the  higher  fluence 
pulses,  containing  0.12  J/cm2  in  each  pulse,  was  slightly 
greater  than  the  ablation  threshold  for  picosecond  pulses  on 
metal.  It  is  believed  that  in  the  experiments  of  March  12th, 
the  laser  was  heating  the  samples  to  the  melting  point, 
whereas  during  the  experiments  of  March  23rd,  each  pulse  was 
ablating  a  small  amount  of  material  immediately.  [Ref.  22] 

1 .  Approximation  Methods 

Two  approximation  methods  were  developed  to  help 
determine  the  required  energy  to  be  delivered  to  the  target. 
Both  methods  assumed  that  the  energy  was  being  delivered  at 
a  high  enough  irradiance  to  overcome  any  thermal  diffusion. 
Results  from  estimates  and  earlier  experiments  tend  to 
indicate  that  an  irradiance  level  of  approximately  10  kW/cm2 
will  be  required  to  melt  through  most  missile  materials. 
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2 .  Scaling 

The  TJNAF  FEL,  using  several  hundred  watts  of 
continuous  average  power  was  used  to  simulate  the  damage 
from  a  MW-class  weapon  by  focusing  the  beam  to  a  smaller 
spot  size  to  increase  irradiance.  The  goal  of  these 
experiments  was  to  begin  to  develop  scaling  rules  that  will 
reliably  predict  the  damage  of  a  larger  laser  without  having 
to  bare  the  enormous  cost  of  building  the  large  laser  first. 
The  experiments  showed  that  the  requirement  for  the  laser 
spot  to  be  larger  than  the  thermal  diffusion  length  is  a 
valid  requirement.  The  burn  through  rate  increases 
dramatically  as  the  laser  beam  radius  becomes  larger  than 
the  thermal  diffusion.  Some  of  the  experiments  indicate 
that  care  must  be  taken  when  calculating  the  thermal 
diffusion  length,  since  this  value  may  change  as  the 
material  changes  states . 

3.  FEL  Pulse  Format 

The  extremely  short,  sub-picosecond  pulse  length  of  the 
FEL  beam  is  a  result  of  the  short  electron  bunches  formed  as 
they  are  accelerated  in  RF  cavities.  The  TJNAF  FEL  has  a 
unique  pulse  format  compared  to  conventional  lasers  with  a 
rapid  sequence  of  short,  powerful  pulses.  The  peak  power  in 
each  pulse  is  about  100  MW  lasting  for  only  about  one-half 
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picosecond  coming  at  a  rate  of  37  MHz.  Other  studies  have 
shown  that  such  ultra-short  pulses  may  give  as  much  as  a 
factor  of  ten  advantage  in  the  reduced  fluence  required  to 
produce  damage  [Ref.  20].  The  experiments  reported  in  this 
chapter  began  to  collect  data  intended  to  show  whether  this 
advantage  exists.  The  results  indicate  that  for  a  fixed 
average  power  level,  a  lower  pulse  repetition  frequency  with 
the  corresponding  higher  power  per  pulse  provides  a  faster 
burn  through  rate.  This  is  probably  because  the  lower  PRF 
pulses  contained  enough  fluence  per  pulse  to  immediately 
vaporize  the  material  as  well  as  raising  it  to  melting 
temperature. 

As  expected  there  is  a  strong  relationship  between  the 
penetration  depth  rate  and  the  laser  irradiance.  When  the 

irradiance  changed  by  a  factor  of  20  (500  W/  cm2  — »  10 
kW  /  cm2 )  the  penetration  depth  rate  changed  by  almost  the 
same  factor. 

From  earlier  studies  and  the  experiments  reported  in 
this  chapter,  it  appears  that  there  is  a  possibility  of 
decreasing  the  fluence  required  to  cause  damage  to  a 
material  with  ultra-short  picosecond  pulses.  This  would  be 
a  most  welcome  finding  for  the  high  energy  laser  community 
since  more  power  means  more  fuel,  more  weight,  and  a  higher 
possibility  of  thermal  blooming  of  the  beam  as  it  propagates 
through  the  atmosphere. 
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4 .  Future  Experiments 

The  TJNAF  FEL  is  scheduled  for  an  upgrade  to  10  kW  of 
power.  This  increase  will  allow  more  flexibility  in  scaling 
experiments  and  further  tests  of  scaling  itself.  Additional 
experiments  that  include  weighing  of  the  samples  before  and 
after  each  run,  new  wavelengths,  changing  wavelengths  during 
irradiation,  new  pulse  formats,  and  other  sample  materials 
would  be  beneficial.  As  experimental  procedures  are  refined 
and  the  amount  of  data  increases,  a  more  thorough  analysis 
of  the  FEL  damage  results  and  comparison  to  other  lasers 
will  become  possible.  Additional  experiments  can  explore 
just  how  high  the  laser  irradiance  needs  to  be  to  bring 
various  missile  materials  to  vaporization  temperature. 

When  higher  laser  power  is  available  in  the  future,  the 
same  irradiations  should  be  conducted  with  bigger  spot  radii 
in  order  to  compare  the  results  and  establish  a  scaling  law. 
It  is  suggested  that  thicker  samples  should  also  be  tested 
in  order  to  determine  if  the  penetration  rate  stays  the  same 
after  the  laser  beam  has  penetrated  the  material  several 
millimeters.  Samples  should  be  machined  to  be  only  slightly 
larger  than  the  laser  spot  so  that  thermal  diffusion  cannot 
occur.  Machining  the  samples  to  resemble  a  small  nose  cone 
shape  may  also  be  important,  because  different  angles  of 
incidence  would  be  included. 


Furthermore,  a  procedure  for  making  the  measurements  of 
burn  through  times  more  accurate  should  be  established,  as 
they  are  very  important  in  determining  penetration  depth 
rates . 
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III.  BASIC  FEL  THEORY 


The  free  electron  laser  (FEL)  has  three  primary 
components:  a  relativistic  electron  beam,  an  undulator,  and 
an  optical  cavity.  The  energy  from  the  electron  beam  is 
coupled  to  the  electromagnetic  field  stored  in  the  optical 
cavity  to  produce  a  very  high  power  coherent  laser  beam.  A 
representation  of  all  three  primary  components  of  an  FEL  is 
shown  in  Figure  5. 


electron  beam  drift  space 


undulator  magnets 


Figure  5.  FEL  Undulator  in  a  Resonator  Cavity 

Figure  5  also  shows  a  drift  space  that  is  included  in 
some  undulator  designs.  The  drift  space  allows  increased 
gain  in  the  second  section  of  undulator  magnets.  This 
design,  known  as  a  klystron  undulator,  will  be  discussed  in 
chapter  IV. 
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One  of  the  most  important  qualities  of  the  FEL  is-  the 
use  of  a  relativistic  beam.  Since  the  electrons  are  moving 
at  speeds  close  to  the  speed  of  light,  they  experience  a 
Lorentz  contracted  undulating  magnetic  field,  which  results 
in  the  emission  of  electromagnetic  pulses  at  short  optical 
or  X-ray  wavelengths.  [Ref.  26]  The  emitted  wavelength  can 
be  calculated  by  solving  for  the  resonance  condition,  and  is 
shown  later  to  be  approximated  by: 


where  A  =  the  emitted  wavelength,  A0  =  the  length  of  one 
undulator  period,  and  y  =  the  Lorentz  factor  of  the 
relativistic  electron  beam.  A  typical  FEL  will  have  Aa  =  5 
cm,  and  y  =  100,  resulting  in  an  emitted  photons  with 
wavelength  A  =  2.5  |im. 

At  TJNAF,  a  350  keV  photo-cathode  injector  generates 
the  electrons.  The  electrons  must  be  accelerated  well  above 
that  energy  to  provide  for  the  short  optical  pulses  desired 
for  a  weapon  class  FEL.  The  electrons  leave  the  injector 
and  proceed  through  a  RF  linear  accelerator,  where  they  are 
accelerated  to  energy  levels  of  40+  MeV. 
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After  exiting  the  accelerators,  the  electron  beam 
enters  the  magnetic  undulator  where  it  interacts  with  both 
the  static  magnetic  field  of  the  undulator,  and  the  dynamic 
optical  field  transiting  through  the  undulator  at  the  speed 
of  light.  The  undulator  can  be  helical  or  linear,  with  the 
helical  undulator  producing  circularly  polarized  light  and 
the  linear  undulator  producing  linear  polarized  light.  The 
electrons  are  moving  at  speeds  very  close  to  the  speed  of 
light,  and  are  almost  in  synchronism  with  the  light  wave. 
However,  special  relativity  requires  that  the  finite  energy 
electrons  travel  slower  than  the  light  beam.  A  condition 
called  resonance  is  defined  as:  the  electrons  are  traveling 
at  a  speed  such  that  over  the  distance  of  one  period  in  the 
undulator,  a  light  wave  of  wavelength  X  will  overtake  the 
electron  beam  by  exactly  one  wavelength.  At  resonance,  the 
effects  of  the  moving  electrons,  the  optical  wave,  and  the 
static  magnetic  field  of  the  undulator  produce  a  beat  wave. 
This  beat  wave  travels  slower  than  the  light  wave  and  is 
called  a  ponderomotive  wave.  At  resonance  the  electrons  are 
actually  moving  at  the  same  speed  as  the  ponderomotive  wave, 
and  the  electric  field  and  resulting  force  the  electrons 
experience  is  constant  and  determined  by  the  portion  of  the 
ponderomotive  wave  that  they  are  traveling  with.  It  is  this 
ponderomotive  force  that  causes  the  electrons  to  slow  down 
or  speed  up.  [Ref.  27] 
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The  change  in  speeds  of  the  electrons  depending  on 
their  phase  with  respect  to  the  ponderomotive  wave  causes 
the  electrons  to  bunch.  The  bunching  is  critical  if  the 
electrons  are  to  add  coherent  radiation  to  the  light  beam. 

The  undulator  is  made  from  either  permanent  magnets  or 
electromagnets  arranged  to  create  an  alternating  magnetic 
field.  The  polarization  can  be  helical  or  linear,  linear 
being  more  common.  The  purpose  of  the  undulator  is  to 
periodically  "wiggle"  the  electrons  as  they  transverse  the 
length  of  the  undulator.  This  wiggling  motion  causes  the 
electrons  to  speed  up  or  slow  down  (absorb  or  emit 
radiation) ,  depending  on  which  part  of  the  ponderomotive 
wave  the  electron  is  in  as  it  is  deflected  by  the  undulator. 

As  with  most  lasers,  the  optical  cavity  stores  the 
light  produced  by  the  laser,  allowing  a  fraction  to  escape 
on  each  pass  of  the  light  wave.  The  amount  of  light  that 
escapes  on  each  pass  is  usually  determined  by  one  of  mirrors 
having  a  slightly  less  than  perfect  reflection  coefficient. 
Allowing  too  much  light  to  escape  results  in  not  enough  gain 
for  the  laser  to  operate,  while  allowing  too  little  light  to 
escape  prevents  the  laser  from  achieving  sufficient  output 
power  to  destroy  a  target. 
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A.  RESONANCE 

Resonance  is  described  above  as  the  condition  where 
over  the  length  of  one  period  of  the  undulator,  one 
wavelength  of  light  passes  over  an  electron.  Over  N  such 
periods  of  the  undulator,  the  light  wave  should  be  a 
distance  of  nA  ahead  of  the  reference  electron.  This 
distance  is  referred  to  as  the  slippage  distance,  since  it 
represents  the  distance  that  an  electron  has  slipped  behind 
the  light  wave.  From  this  information  the  wavelength  of  the 
light  can  be  calculated.  First,  the  total  distance  the 
light  wave  moves  ahead  of  the  electrons  is  given  by: 

N/l  =  (c  —  v2)  At 

N/l  =  (c  —  c&)  At  (21.) 

N>3.  =  c  (l  -  Pz)  At 

where  N  is  the  number  of  undulator  periods,  A  is  the 
wavelength  of  the  light  wave,  c  =  the  speed  of  light,  vz  = 
the  speed  of  the  electrons  along  the  line  of  the  undulator, 
Pz  =  vz/c,  and  At  =  N^v*  =  the  time  for  the  electrons 

to  transverse  N  undulator  periods.  Simplifying  equation 
(21.)  to  solve  for  the  wavelength  in  terms  of  velocity  in 
the  longitudinal  direction: 
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(22.) 


NA  =  c  (l  -  A) 


NA0 

Ac 


A  = 


i- A 

A 


An  approximation  of  the  wavelength  can  be  made  by  using 
a  binomial  expansion  to  determine  /?z.  The  transverse 

K 

velocity  is  shown  later  to  be  f$L  =  — - ,  where  Ku  =  eBAo/2;zmc2 

Y 

is  the  dimensionless  undulator  parameter.  Included  in  Ku 
are  e,  the  electron  charge  magnitude,  B,  the  undulator  root- 
mean-square  field  strength,  Ac  the  undulator  wavelength,  m 
the  mass  of  the  electron,  and  c  the  speed  of  light.  Ku  is  a 
dimensionless  number  that  gives  a  general  measure  of  the 
magnetic  field  strength.  Continuing  with  the  derivation  for 
optical  wavelength  with  the  velocity  of  the  electrons  with 

respect  to  the  velocity  of  light,  (3, 


r  = 


i 


1  -  l 3 2 

r 2  =  i 


i  -  A2  -  A2 


-  A 


(23.) 

(24.) 


r_2(i  +  K2)  =  1  -  pi 


(25.) 
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PI  =  1  - 


(26.) 


fe  +  *u) 

r2 


A 


1  +  K. 

2r2 


2  '\ 


(27.) 


where  the  last  approximation  was  made  by  using  the  binomial 
expansion  since  Ku  is  typically  of  order  unity,  and  y2  of 
order  104.  Inserting  this  approximation  for  fiz  back  into 
equation  (22.),  the  value  for  the  optical  wavelength 
becomes : 


X  = 


1  -  A 
A 


A  = 


( 

1  -  1  + 

'i  +  KT\ 

_ 

2y2  Jj 

/ 

1  +  KuT) 

1  - 

V 

.  2 y2  JJ 

(28.) 


The  term  in  square  brackets  is  very  small  and  can  be 
dropped  from  the  denominator.  However,  it  is  the  only 
remaining  term  in  the  numerator  and  must  be  kept. 
Therefore,  the  optical  wavelength  becomes: 


45 


1  +  K* 


(29.) 


/I  = 


2  y 


and  since  K  is  usually  <1, 


A  * 


2  y 


(30.) 


Equation  (29.)  shows  how  the  optical  wavelength  can  be 
controlled  by  three  parameters,  the  undulator  wavelength, 
the  undulator  parameter,  which  is  primarily  controlled  by 
the  strength  of  the  magnetic  field,  and  the  energy  of  the 
electron  beam. 

The  resonance  condition  used  to  calculate  the 
approximate  wavelength  of  the  light  beam  is  actually  not  the 
desired  state  for  the  electron  beam.  At  resonance  just  as 
many  electrons  would  be  absorbing  energy  from  the  optical 
beam  and  giving  energy  back  to  the  optical  beam.  At  low 
gain,  the  optimum  gain  condition  is  for  the  average  electron 
phase  velocity,  vo,  with  respect  to  the  ponderomotive  wave 
to  be  at  vo  =  2.6,  which  means  the  electrons  are  initially 
traveling  slightly  faster  than  the  ponderomotive  wave.  This 
condition  will  cause  the  electrons  to  bunch  on  the  positive 
slope  of  the  ponderomotive  potential,  which  will  cause  them 
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all  to  slow  down  (give  up  energy)  together  resulting  in 
coherent  photon  emission. 

B.  ELECTRON  MOTION 

Three  fields  combine  to  form  the  ponderomotive  wave, 
producing  three  forces  on  the  electrons  as  they  move  through 
the  undulator.  Two  of  the  forces  are  magnetic.  They  are 
caused  by  the  alternating  static  magnetic  field  of  the 
undulator,  and  the  moving  magnetic  field  of  the  optical 
wave.  The  third  force  is  from  the  moving  electric  field  of 
the  optical  wave.  For  a  helical  undulator  the  three  fields 
written  in  CGS  units  are: 

Bu  =  Bu  (cos(k0z),  sin(k0z),  0)  (31.) 

Es  =  Es  (cos(^),  -  sin(^),  0)  (32.) 

and 

Bs  =  Es  (sin((^),  cos(^),  0)  (33.) 

where  Bu  is  the  magnetic  field  of  the  undulator,  Es  and  Bs 
are  the  electric  and  magnetic  portions  of  the  optical  field, 
k  is  the  wave  number  of  the  undulator,  and  l//=  (kz-fyt+0)  is 

o 

the  phase  of  the  optical  wave,  where  <|>  is  the  initial  phase 
at  t  =  0  and  z  =  0 . 
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The  total  force  acting  on  the  electron  is  given  by  the 
Lorentz  force  equations  [Ref.  28,  pg.  579].  The  equations 
are: 


=  -  —  (is  +  P  x  i)  (34 . ) 

dt  me 

and 

^  is  (35.) 

dt  me 

Equations  (34.)  and  (35.)  are  actually  four  equations. 
Using  these  four  equations  plus  equation  (23.) 


r2  = 


i  -  ft  i- ft-  ft 


(23.  ) 


gives  a  total  of  five  equations  and  four  unknowns.  As  is 
commonly  done,  since  the  motion  of  interest  is  in  the 


transverse  direction,  the 


dQffz) 


dt 


equation  will  be  ignored. 


The  transverse  components  can  be  written  separately  as 
components  in  x  and  y  as: 
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^  =  -(h„*Axb„  +  AxbJ 

dt  me 

=  — -  (eb  cos(^)  -  Bu/?z  sin(k0z)  -  Es/?z  cos(l//j)  (36.  ) 

dt  me 

d-—  =  —  (Es(l  -  fiz)  cos(^)  -  Bu/?z  sin(kcz)) 
dt  me 


and 


=  —  (Esy  +  Pz  x  Bux  +  £  x  Bsx ) 
dt  me 

=  — -  (-  Es  sin(yr)  +  Bu/?z  cos(k0z)  +  E S/3Z  sin (^))  (37 . ) 

dt  me 

d^y -  =  —  (~  Es(l  -  0Z)  sin(yr)  +  Bu/?z  cos(k0z)) 
dt  me 


Since  l-/3z«l,  equations  (36.)  and  (37.)  can  be 
simplified  to: 


d(/& )  =  _ £(_  Bu/?z  sin(k0z))  (38.) 

dt  me 

and 


dt 


— (Bu&cos(k0z)) 
me 


which  can  be  written  together  as 


(39.) 
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(-  sin(koz),cos(koz),0) 


(40.) 


d(r/U_~eB  A 

dt  me 


where  the  zero  has  been  added  to  maintain  the  longitudinal 
dimension  coordinate.  By  inspection  equation  (40.)  can  be 
integrated  to  show  that : 


PL  =  — —  (cos(k0z),  sin(koz),0) 


(41.) 


The  second  Lorentz  equation  is: 


dy  _  —  e 
dt  me 


P  ■  Es 


(42.  ) 


and  using  equation  (32.),  this  can  be  written  as 


d  y 
dt 
dy 
dt 


-  e 


me 


(fiy  fiy  A)  •  (E.  cos  iff-  Es  sin  ^,0) 


-  e  (  -  K„ 


me 


Es  (cos(k0z)  cos  -  sin(kcz)  sin  y. /) 


(43.) 


Now  using  the  trigonometric  identity 
cos (a+b) =cos (a) cos (b) -sin (a) sin (b) ,  this  equation  simplifies 
to 
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(44.) 


—  -  eKljEs  cos(kQz  +  y/) 
dt  me)' 

* Y-  =  cos(C  +  *) 

dt  me)' 


where  £  =  (k  +  k0)z  -  cot,  the  electron  phase  in  the 
ponderomotive  wave. 

Next,  the  time  derivative  of  equation  (23.)  can  be  used 
to  get  a  relation  between  the  time  derivative  of  electron 

y  t  and  the  time  derivative  of  the  fractional 
velocity,  /?  , 

r2{ i  +  *u)  =  i  -  Pi 

-  2r3  r  (i  +  Ku)  =  -2A  A  <45-) 

r  =  r2A  A 

y  (i  +  k*) 

With  (0  =  kc,  the  electron  phase  in  the  ponderomotive 
field  and  its  derivatives  are: 


C  =  (k  +  k0)z(t)  "  ^ 

C  =  (k  +  kc)z  -  0) 

i  =  (k  +  kj&c  -  69 

C  =  (k  +  kj  A  c 
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Therefore, 


z 


•  • 

£ 

(k  +  kQ)c 


(47.  ) 


Now,  inserting  equation  (47.)  into  equation  (45.), 


t  =  y%  k 

r  (i  +  K) 

r  =  r2&  V 

Y  (l  +  Ky)  (k  +  kj  c 


(48.) 


and  since  Pz~l,  k  »  k0 ,  and  kc  =  Q) ,  this  becomes 


•  •  • 

r  _  f  C 

Y  (l  +  <>y 


(49. ) 


-  ,  _  2  7TC  .  . 

Now  using  equation  (29.)  and  replacing  A  with  -  and  Aq 

CO 


with 


2no 

a 


CO  = 


1  +  K* 


(50.) 


so  that  equation  (49.)  becomes: 
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Finally,  using  equation  (44.)  for  y  : 


•  • 

y  =  — —  y  =  eKuEs  cos((  +  0 ) 

2^0  me  y  (52.) 

V*  =  2«0  ^4  cos(C  +  0) 

m  cy 


Equation  (52.)  gives  the  electron's  phase  dynamics  in 
the  form  of  a  pendulum  equation.  It  governs  the  electron's 
phase  space  motion  resulting  from  the  ponderomotive  wave. 
Whenever  the  total  phase  ( £  +  0)  of  an  electron  falls 

between  ±  ,  the  cosine  term  will  be  positive,  resulting 
in  an  acceleration  of  that  electron.  When  the  phase  is 
between  ^  and  ^n/<  the  acceleration  is  negative.  This 

positive  and  negative  acceleration  results  in  the  bunching 
of  the  electrons  on  a  scale  size  of  the  optical  wavelength. 
The  electrons  longitudinal  motion  is  approximately  limited 
to  ±A  .  Introducing  a  the  dimensionless  variable 


__  ct 
L 

where  c  is  the  speed  of  light,  t  =  time,  and  L  is  the  length 
of  the  undulator.  The  pendulum  equation  can  be  rewritten  as 
a  second  order  derivative  with  respect  to  T  as 


^  =  |a|  cos((  +  <j>)  (53  . ) 

where  dr  =  -  dt  and  lal  =  4^Ku^Es  .  a  =  J^e1*  is  called  the 
c  romc 

dimensionless  optical  field.  From  equation  (53.)  it  can  be 
seen  that  the  electron  motion  is  a  function  of  the  optical 
field. 

C.  OPTICAL  WAVE  EQUATION 

In  the  absence  of  the  optical  field,  the  electrons 
experience  the  forces  from  the  undulator  field  alone.  This 
field  causes  the  electrons  to  undergo  small  transverse 
accelerations,  resulting  in  spontaneous  transmissions  of 
photons  in  the  forward  direction.  The  number  of  electrons 
in  a  single  pulse  is  on  the  order  of  109,  and  the  number  of 
photons  emitted  during  the  first  pass  through  the  undulator 
is  enough  that  a  classical  field  is  formed  and  then 
amplified  though  stimulated  emission.  [REF.  29,  p.  129] 


54 


The  optical  wave  equation  written  in  terms  of  the 


vector  potential  is:  [Ref.  28,  pg.  240] 


(54.) 


where  Ji  is  the  transverse  current  density. 

The  vector  potential  for  a  circularly  polarized  plane 
wave  can  be  written  as  a  function  of  z  and  t  only: 

a(z,  t)  =  E^-Z' [cos(t4  -  sin(<4  0]  (55.) 

k 


Carrying  out  the  spatial  and  time  derivatives  of 
equation  (55.),  one  obtains 


d2A  _  1  3es 
dz2  k  dz 


k  + 


V 


[cos(f),-  sin(^),0]  +  [sin(^),  cos(^),0] 

dz  k  dz 


+  -  fk  +  [cos(f),-  sin(^),0]  +  ^  [cos(^),~  sin(^),0] 

k  dz  ~ 


V 

E  (  dth  \2 


k-3 

dz 


dz J  l"_vw'  — —  k  dz2 

[-  sin(^),-  cos(^),0] 


(56.) 


and 
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O)  d2A  dE _ 


at2 

dEs 

at 


at 


(O 


d<(> 

at 


a2E, 


[co s{y/)-  sin(^),0]  +  — r  [sin(^),  cos(^),0] 


at2 


-  CO 


J 


d<t> 

at 

s 

d2(f) 

|[cos(^)  -  sin(^),0]  +  Es  [cos(^),-  sin(^),0] 


V 

+  Esfft  ~  sin(^),_  cos(^)'0] 


(57.) 


Assuming  that  the  optical  phases  and  amplitudes  vary 
slowly  in  time  and  space 


«  kEs,  «  k0,  «  (OEs,  and  «  co<t> 

dz  oz  ut  at 


the  second  order  derivatives  are  very  small  and  can  be 
eliminated.  This  allows  equation  (54.)  to  be  written  as: 


— -2  1  d2  ' 

-h  \ 

ra es  i  aEs^| 

V 

c2  at2_ 

A\x,  t)  =  2 

az  c  at  j 

+  2E_ 


[cos(f),  sin(^),  0] 


—  +  —  -^-1  [-  sin(^),  cos(^),  0]  =  -  —  Ji (x,  t) 

az  c  at  J  c 


(58.) 


At  this  point  it  is  convenient  to  use  the  vectors  of 
equation  (58.)  to  define  two  unit  vectors: 


s  1  =  [cos(^),  sin(^),0]  and  £i  =  [-  sin(^),  cos(^),0] 
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The  wave  equation  (58.)  can  now  be  separated  into  two 
orthogonal  wave  equations : 


and 


rdEs  +  J_  dEs' 
^dz  c  3t  J 


(59.) 


d0 

dz 


1 

c  dtj 


2n 

c 


(j  i 


A  > 
£  2 

) 


(60.) 


Still  left  to  be  determined  is  the  transverse  current. 
It  can  be  written  as  the  sum  of  individual  electron 
currents : 


J±(x)  =  -|e|c^  fi±S3(x.  -  ri)  (61.) 

i 


where  r±  is  the  position  of  the  ith  electron  with  charge  e 

and  transverse  velocity  ft Lc ,  and  8  is  the  three  dimensional 

Dirac  delta  function.  Substituting  equation  (41.)  into 
equation  (61.)  gives: 


J±\x 


=  -|e|  c£  — —  (cos(k0z),  sin(k0z),  0)  £3(x  -  ri)  (62  . ) 
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which  can  be  projected  onto  the  two  unit  vectors  defined 
above  to  obtain: 


3x0.  a,  =  14*. s co<¥  +  kA  -  n)  <63.) 

i  Y 

3x0  •  i,  =  |eK2  Sin(y  +  ^  -  n)  <64. ) 

i  Y 

When  the  number  of  undulator  periods  is  very  large,  N 
»1,  then  the  gain  bandwidth  curve  is  very  narrow.  The 

change  in  energy  of  the  electrons  per  pass  through  the 

undulator  is  very  small  and  therefore  the  Lorentz  factor  can 
be  pulled  out  of  the  integral . 


.  Ji(x)  •  £i  =  ^CKu-  ^  cos(^  +  kQz)  <53(x  -  ri)  (65.  ) 

Y  i 

Jx(x)  •  £ 2  =  ^CKu  sin (y/  +  kQz)  J3(x  -  n)  (66. ) 

Y  i 

The  sum  of  the  product  of  the  individual  electrons 
times  the  phase  of  that  electron  in  the  ponderomotive 
potential  can  be  replaced  by  considering  a  volume  element  of 
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electrons  of  density  pe  times  the  average  phase  of  the 
potential  in  that  volume  element. 


e|cKupe 

r 


(cos{ys  +  kQz)) 


(67.) 


and 


McKupe 

r 


(sia(w  +  kQz)) 


(68.) 


Now,  rewriting  y/  +  k0z  as  C  +  0  >  the  two  orthogonal 
wave  equations  (59.)  and  (60.)  become: 


and 


'dE,  +  I 

,  dz  c  9t  , 


2fr|e|Kupe 

r 


(cos (C  +  0)) 


(69.) 


,  (d0_  t  1 

s^8z  c  dt 


2#  upe 

r 


(sin(C  +  *)) 


(70.  ) 


Equations  (69.)  and  (70.)  can  be  formed  from  the  real 
and  imaginary  parts  of  the  following  equation: 


d  ,  1  3  V  2*Hk» ft 

_ —  —  £j  e  — 

3z  c  3t  J  / 


<e-> 


(71.  ) 
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Assuming  a  long  electron  pulse,  there  will  be  no 
spatial  dependence  in  the  z  direction,  and  no  longitudinal 

d 

optical  modes  are  followed.  Therefore  the  —  term  can  be 

dz 

dropped.  [Ref.  29,  p  132] 


I  a  (E,e,)  .  - 

c  dt  Y 


(72.) 


Finally,  using  the  dimensionless  variables 


a  = 


_  4;r|e|KuLEse 


i0 


Y  me 


and 


8N(e^KuL)2pe 

/3mc2 


the  optical  wave  equation  can  be  written  as 


(73.) 


where  a  is  the  derivative  of  the  dimensionless  optical 
field,  a,  with  respect  to  the  dimensionless  time  variable  x. 
Equation  (73.)  shows  that  with  no  current,  or  with  zero 
average  phase  change  (no  blanching)  ,  there  is  no  change  in 
the  optical  field. 
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D.  SHORT  PULSES 

Recall  that  at  resonance,  the  slippage  distance  was 
defined  as  nA.  Short  pulses  are  defined  as  pulses  whose 

length  is  less  than  or  equal  to  the  slippage  distance. 

In  the  previous  section,  the  long  pulse  assumptions 
were  used,  basically  stating  that  because  the  pulse  length 
was  larger  than  the  slippage  distance,  the  optical  pulse 
continued  to  see  the  same  electron  density  throughout  the 
transition  through  the  undulator  over  many  passes.  For 
short  pulses  this  assumption  cannot  be  made. 

Two  processes  are  occurring  simultaneously  during  the 
passage  through  the  optical  cavity.  The  first  process  is 
the  electrical  pulse  slipping  backwards  with  respect  to  the 
optical  pulse.  The  second  is  the  bunching  that  is  occurring 
on  the  optical  wavelength  scale.  Due  to  the  slippage  of  the 
short  electron  pulses,  only  the  rear  section  of  the  optical 
pulse  will  remain  overlapped  with  the  electron  pulse  by  the 
time  the  electrons  have  achieved  bunching.  This  results  in 
the  optical  pulse  amplification  taking  place  only  in  the 
rear  of  the  pulse,  while  the  entire  pulse  experiences  the 
losses  associated  with  the  cavity.  The  pendulum  equation 
showed  that  electron  motion  (phase  change)  only  takes  place 
in  the  presence  of  the  optical  field  and  is  proportional  to 
the  magnitude  of  the  optical  field. 
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a  cos(^  +  <[>) 


(53.) 


(T 


Since  the  amplification  is  only  taking  place  on  the 
rear  portion  of  the  optical  pulse,  the  centroid  of  the  pulse 
grows  toward  the  rear  (making  the  light  pulse  appear  to  be 
moving  slower  than  c)  with  each  pass  through  the  undulator. 
Eventually,  the  slower  moving  optical  pulse  is  entirely  off 
the  electron  pulse,  and  no  bunching  takes  place.  The 
optical  pulse  eventually  decays  due  to  the  cavity  losses. 
This  effect  is  shown  in  the  left  column  of  Figure  6.  The 
center  window  of  this  column  shows  the  optical  field, 
a ( z , n)  ,  growing  strong  at  n  =  100,  and  then  decaying  again 
around  n  =  350. 


-1.5  z  1.5-47  V  47  -47  V _ 47 


-1.5  z  1.5-47  V  47  0  n  400 


Figure  6.  Simulation  of  Short  Pulse  Evolution  in  the  FEL. 
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Figure  6  shows  9  sections  arranged  in  a  3  x  3  grid 
pattern.  The  lower  left  corner  indicates  the  position  of 
the  parabolic  electron  pulse  shape.  The  electron  pulse 
slips  back  by  the  slippage  distance  NX  relative  to  the  light 
pulse  as  both  pulses  move  through  the  undulator  from  T=  0 
to  x  =1 .  The  middle-left  plot  shows  the  evolution  of  the 
optical  field  |a(z,  n)|  as  it  grows  over  n  =  400  passes  through 
the  undulator,  with  the  final  optical  pulse  shape  shown  in 
the  top- left  window. 

The  center  graph  shows  the  evolution  of  the  optical 
power  spectrum,  P(v,n)  over  n  =  400  passes  through  the 
undulator,  where  v  is  the  electron  velocity  or  energy.  The 
top-center  plot  shows  the  final  optical  power  spectrum  after 
n  =  400  passes.  The  bottom-middle  plot  is  the  weak-field 
gain  spectrum,  G(v) ,  showing  a  maximum  gain  occurring  at  v  = 
2.6. 

The  lower-right  window  shows  the  Power  of  the  optical 

pulse  P(n).  The  middle-right  plot  shows  the  evolution  of 

« 

the  electron  spectrum  f  (v,n)  over  n  =  400  passes,  while  the 
top-right  plot  shows  the  final  spectrum  after  n  =  400  passes 
through  the  undulator. 

The  electron  pulse  is  timed  to  enter  the  cavity  in 
synchronism  with  the  reflected  optical  pulse,  which  has  a 
time  interval  of  2S/c,  where  S  is  the  one  way  distance 
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between  the  mirrors  of  the  cavity.  However,  as  the  optical 
pulse  starts  to  move  back,  the  bunching  of  the  electron 
pulse  is  delayed  more  and  more,  causing  the  amplification 
process  to  take  place  even  further  back  on  the  optical 
pulse.  Eventually,  the  overlap  between  optical  and  electron 
pulses  is  zero,  and  the  optical  pulse  decays  to  zero. 

The  desynchronism,  d,  is  defined  as  the  displacement 
between  the  electron  and  optical  pulses  at  the  entrance  to 
the  optical  cavity.  It  is  normalized  to  the  slippage 
distance.  d  =  0  means  the  pulses  are  in  perfect 
synchronization.  However,  since  in  effect  the  optical  pulse 
is  moving  back  with  each  pass,  the  desynchronism  needs  to  be 
adjusted  so  that  the  "slowed"  optical  pulse  arrives  with  the 
electron  pulse.  This  is  done  by  shortening  the  length  of 
the  optical  cavity  by  an  amount  AS,  resulting  in  a 

desynchronism  value  of  d  =  2AS/NA. 


E.  GAIN 

FELs  are  generally  operated  as  either  amplifiers  or 
oscillators.  Amplifiers  are  single  pass  FELs  that  require 
high  gain  to  develop  high  power  in  one  pass  through  the 
undulator.  Oscillators  create  their  power  over  many  passes 
and  therefore  can  operate  at  low  gain. 
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The  gain  of  a  FEL  is  defined  as  the  percentage  change 
in  the  optical  field  per  pass  through  the  oscillator: 


G  = 


(74.) 


where  ax  is  the  optical  field  strength  at  the  end  of  the 
undulator,  and  a0  is  the  optical  field  strength  at  the 
beginning  of  the  undulator. 

To  achieve  gain,  energy  has  to  be  removed  from  the 
electron  beam  and  added  to  the  optical  beam.  Therefore  a 
basic  method  to  analyze  the  gain  in  the  optical  field  energy 


is  to  determine  the  change  in  energy  of  the  electrons. 


Ay 

r 


for  each  electron,  and  add  these  changes  together  for  all 
the  electrons.  A  negative  change  in  average  energy  of  the 
electron  beam  means  a  positive  change  in  the  energy  of  the 
optical  beam. 

The  change  in  energy  of  a  single  electron  is 
proportional  to  its  change  in  phase  velocity.  The 
instantaneous  phase  velocity  is  given  in  equation  (46.)  as: 


v 


=  i  =  (k  +  k)Pz°  -  <° 


(46.) 


or  in  units  of  dimensionless  time  z  where  dt  — >  —  dr 

L 
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(75.) 


v  =  C  =  -  I(k  +  kj&c  -  0)} 

(k  +  k0  )/?2 


V  =  C  =  L 


v  =  ^  =  L[(k  +  kj/?z  -  k] 


where  L  =  the  length  of  the  undulator. 

The  change  in  the  phase  velocity  can  be  written  as: 


Av  =  A  C  =  L(k  +  kjA/?z 
Av  =  A  l  ~  LkA J3Z 


where  the  approximation  comes  from  the  fact  that  k  »  ko. 

To  relate  Av  to  Ay  use  the  substitutions  L  =  N/lo,  k  = 

2n/X,  kc  =  2ji!Xo,  and  equation  (29.) 


X  = 


1  +  K* 

2  y2 


(29.) 


to  obtain: 
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Av 

Av 

Av 


n2, 

N2tt| 


f  2 7t^ 


°U; 

(X\ 
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N2/r 


2  Y 


m 

2  A 


1  +  K2 


m 


(77.) 


Now  use  equation  (48.)  in  difference  rather  than 
differential  form: 


A  y  _  fA  Ayg 
r  (i  +  Ku) 


to  substitute  A/?z  into  equation  (77.) 


Av 

Av 


N2d 


2/2 


1  +  K2 


4^N  A  y 


A Y  &  +  K2) 

r3  A 


A  Y 


Av 


Am 


Ay 


r 


(78.) 


where  the  last  approximation  comes  from  the  fact  that 
Pz  =  1 .  The  change  in  energy  of  the  electron  is  Ajmc2or 
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Av  2 

AEC  =  - me 

s  4?W 


(79.) 


The  average  energy  change,  AEs,  for  an  electron  is: 


AES  =  A/Totalmc 


(Av)  2 
- — -  me 
4;zN 


(80.) 


where  A/Total  is  the  Siam  total  change  in  energy  of  the 

electron  beam  from  all  the  electrons. 

Assuming  a  monoenergetic  beam  enters  the  undulator  with 
phase  velocity  vo ,  this  can  be  rewritten  as: 


.  (v)  -  vo  , 

AEs  =  A/Totalmc2  =  —  /me 


(81.) 


The  final  term  to  determine  is  the  average  energy  term 
for  each  electron  at  the  end  of  the  undulator,  (v)  .  The 

phase  velocity  v  can  be  determined  by  integrating  the 
pendulum  equation  once 


=  j  C  dr  =  J  |a|  cos(C  +  (/>)  dr 


(82.) 
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In  weak  fields  the  changes  in  the  optical  field, 

o 

electron  phase,  and  optical  phase  are  small:  a  =  0  and 

<j>  «  0  .  In  weak  fields  the  pendulum  equation  can  be 
rewritten  by  expanding  the  phase,  £,  in  terms  of  the 
dimensionless  optical  field  term  a 

'£  =  -il  (c(0)  +  c{1)  +  C{2)  +  ...)  =  |a|  cos {c{0)  +  c(1)  +  C2)  +  .  •  •)  (83. ) 
dr 

using  the  zeroth  order  conditions 

a(0)  =  a 

o 

v(0)  =  vD  +  .  .  .  (84.) 

C(0)  =  Co  +Vor  +  ... 

To  zeroth  order  the  pendulum  equation  becomes 

O°(0)  .  V 

(  =^r(fw’)  =  o  <85-> 

To  first  order  in  'a'  the  equations  becomes 

C  =  (C(0)  +  C0))  =  |a(0)|  cos(^(0))  =  a0  cos(£o  +  v0t)  (86.) 
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which  from  equation  (85.)  is 

C  =  (c{l})  =  ao  cos(^0  +  V0t)  (87.) 

Integrating  equation  (87.)  once  produces 

C  =  ^(I)  =  —  [sinfe  +  Vor)  ~  sin(Co)]  <88-) 

Integrating  a  second  time  provides  the  first  order 
approximation  for 


C(1)  =  ^j[-  cos{(0  +  v0t)  +  cos(Co)]  -  —  [sin(Co)] 

V°  V°  (89.) 

Cil)  =  ~  ^  !c°s(^0  +  v0T)  ~  cos(^o)]  +  V0t[sin(£,)]} 

Using  these  equations  the  second  order  approximation  of 


the  pendulum  equation  is: 


..(2) 

c 


.0(2) 

c 


00(2) 

c 


■fj  (f  (.)  +  f<»  +  fB>)  =  |a|  cosfe- <">  +  f  ">) 

dr2 


a2  < 


dr 


Co  +  voT  +  -  ~T  {[cos(Co  +  v)  ”  c°s(^o)]  +  v0t[sin(^0)]}  +  C 


(2) 


Ko 
=  lal 


cos(C(0)  +  C") 


^  +  a0  cos(£,  +  v0r)  +  {C{2))  =  a0  cos(C0  +  vj  +  C(0) 


(90.) 


The  right  side  can  be  expanded  by  the  identity 
cos (a+b) =cos (a) cos (b) -sin (a) sin (b)  and  by  the  approximations 

cos(^(1))  «  1  and  sin(^(1))  =  C{1)  '• 


a0  cos{C0  +  v0t)  +  (C{2))  =  a0[cos(^0  +  v0r)  cos(C{l))  -  sin(£,  +  v0t)  sin 

OT 

a0  cos(£o  +  voT)  +  TT  (C{2))  =  ao[cos(£o  +  voT)  ~  sin(Co  +  vOfe'0*)] 

OT 

■jT  (C{2))  =  _ao  sin(^0  +  v0r)(^(1)) 

OT 

TT  (C{2))  =  %  sin(£0  +  v0rXcos(C0  +  v0T)  ~  cosiCo )  +  vofc  sin(Co)} 
dr  v0 

(91.) 


Integrating  equation  (91.)  once  gives: 


v(2)  =  £2. 
^0 


1  -  V0r  sin  Co  cos (Co 


(92.  ) 
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Equations  (89.)  and  (93.)  give  higher  order 
approximations  for  the  phase  and  phase  velocity  of  the 
electrons,  which  can  be  represented  together  in  phase  space 
as  shown  in  the  large  window  of  Figure  7 .  This  figure  shows 
the  phase  space  position  of  20  electrons  after  transversing 
the  distance  of  the  undulator.  About  half  of  the  electrons 
have  bunched  together  at  a  phase  position  11,=  n. 


-Jt/2  5  3JI/2  OX  1 


Figure  7.  Simulation  of  Phase  Space  Plot  of  20  Electrons. 

In  Figure  7  the  large  window  shows  the  electrons 
position,  and  velocity,  v,  relative  to  the  ponderomotive 

wave.  This  example  shows  only  the  final  phase  space 
positions  of  the  electrons.  The  lower  right  window  shows 
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the  evolution  in  the  optical  phase,  <t>(t),  and  the  upper  right 
window  shows  the  gain  G(T)  . 

Equation  (89.)  shows  that  to  first  order  the  bunching 
is  proportional  to  the  optical  field  a0 ,  and  that  some  of 
the  electrons  move  forward  while  some  move  backwards, 
depending  upon  their  initial  phase  position,  £0  and  phase 

velocity,  v0 . 

To  first  order,  the  final  velocity  is  proportional  to 
a0 .  It  is  also  dependent  on  initial  phase  and  phase 
velocity,  with  some  of  the  electrons  gaining  energy  and  some 
losing  energy. 

To  complete  equation  (81.),  the  final  unknown  (v)  can 
be  approximated  as: 


(v)  =  (v(0)  +  v(1)  +  v(2)  +  .  .  .)  (93.) 

The  electrons  enter  the  undulator  with  a  pulse  length 
that  is  many  optical  wavelengths  long.  The  phase  space  is 
periodic  on  the  optical  wavelength  scale,  so  the  average 

change  in  energy  AES  can  be  determined  over  one  wavelength: 

55.  =  <94-> 

i  0  ^ 
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The  average  for  the  zeroth  order  approximation  is: 


(v°)  =  v0  (95.) 

indicating  no  change  in  average  energy  of  an  electron. 

To  first  order  the  average  is: 


(v(1))  =  v0  +  0 


(96.) 


since  the  two  sine  functions  in  equation  (88.)  cancel  when 
averaged  over  all  values  of 

To  second  order  the  average  is: 


.2  r 


vn  +  0  + 


0  L 


cos(v0r)  -  1  +  ^  v0z  sin(v0r) 


(97.) 


which  produces  a  net  change  in  the  average  energy  dependent 
upon  the  initial  phase  velocity  of  the  electron. 

Equation  (97.)  can  be  placed  into  equation  (81.)  to 
give  the  average  phase  velocity  change  in  the  electron  beam. 

The  gain  in  the  optical  beam  was  defined  as  the  change 
in  energy  of  the  optical  beam  divided  by  the  initial  energy . 
It  can  be  determined  on  the  differential  volume  element  dV 
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level  as  the  change  in  the  energy  of  the  electron  beam  in 
that  volume,  divided  by  the  initial  optical  beam  energy  in 
the  same  volume.  The  energy  in  a  volume  element  containing 
the  optical  beam  is: 


-'O-beam 


2e1  . 
— dv 

8 71 


(98.) 


where  E  is  the  Electric  field  of  the  beam. 

The  average  change  in  energy  per  electron  is  given  in 
equation  (97.).  The  number  of  electrons  in  a  volume  element 
of  light  dV  is  pFdV,  where  F  is  the  filling  factor,  which  is 
the  ratio  of  the  electron  beam  area  to  the  optical  beam 
area.  The  gain  in  the  optical  beam  is  therefore: 

(pFdv)[ymc2((v)  -  Vp)  /  47Zn]  (99  } 

2E2dV  /  S;r 

and  using  equation  (97.)  the  gain  becomes: 


G  = 


de. 


e-beam 


^  c  o-beam  J 


G  =  - 


47zpeFymc2  a^ 


cos(v0r)  -  1  +  -  v0r  sin(v0r) 


(100. ) 


Using  the  dimensionless  variable  j  =  8N(e^K2L)2  peF  /  y3mc2 
produces 
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G 


(101. ) 


j[2  -  2  cos(v0t)  -  v0t  sin(v0T)] 


Equation  (101.)  shows  that  in  the  weak  field  low 
current  approximations  used  here,  the  gain  is  primarily  a 
function  of  initial  phase  velocity,  v0 .  Figure  8  below  is  a 

plot  of  the  gain  spectrum  for  a  range  of  phase  velocities 
from  -  15  to  +  15.  Gain  G  is  plotted  on  the  vertical  axis, 
while  initial  phase  velocity  v0  is  plotted  on  the  horizontal 

axis.  The  computer  simulated  spectrum  shows  a  peak  gain  of 
27%  when  the  phase  velocity  equals  2.6.  The  spectrum  is 
anti-symmetric  about  V0  =  0  ,  indicating  that  the  maximum 
loss  would  occur  for  a  electron  beam  with  a  velocity  of 
v0  =  -2.6. 

0.2691 

0,0 

-0.2691 

Figure  8.  Weak  Field  Gain  Spectrum 
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A  phase  space  plot  showing  the  evolutions  in  phase 
space  of  20  electrons  with  initial  phase  velocity  of 
V0  =  2.6  is  shown  in  Figure  9.  The  light  colored  dots 

indicate  the  phase  space  positions  at  the  beginning  of  the 
undulator,  and  the  darker  colored  dots  show  the  positions  at 
the  end.  As  can  be  seen  in  the  figure,  many  of  the  darker 
colored  dots  are  bunched  in  the  phase  position  of 
approximately  C,  —  n ,  and  are  all  moving  down  in  phase 
velocity.  This  indicates  that  by  the  time  the  electrons 
reach  the  end  of  the  undulator,  they  are  bunched  together  on 
the  optical  wavelength  scale,  and  are  losing  energy  to  the 
light  beam. 


Gain  0.3 

Li 


0.1 
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Figure  9.  Phase  Space  Evolution  for  20  Electrons. 
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F.  TAPERED  UNDULATOR 

The  resonant  condition  was  defined  in  equation  (29.) 

as : 


A  = 


1  +  K* 

2  y2 


/L 


(29.) 


The  resonant  condition  was  used  to  calculate  the  change 
in  energy  expression  in  equation  (78.)  .  However,  once  the 
optical  beam  amplitude  becomes  large,  the  change  in  energy 
of  the  electrons  is  such  that  they  are  no  longer  in 
resonance  with  the  optical  beam  and  the  gain  decreases. 

One  method  to  maintain  resonance  and  as  a  result 
maintain  gain  is  to  change  one  of  the  parameters  in  equation 
(29.),  so  that  the  right  side  of  the  equation  remains 
constant  even  with  decreasing  energy.  This  can  be  done  by 
decreasing  the  undulator  wavelength,  Aq,  or  by  decreasing 
the  magnetic  field  strength  resulting  in  a  decreased  value 
in  the  undulator  parameter  Ku.  Both  methods  are  referred  to 
as  tapering  the  undulator. 

Effectively  either  one  of  the  changes  produces  the  same 
effect  on  the  electrons  as  if  the  electron  had  been 
accelerated  back  up  to  the  resonance  speed  by  some 
electrostatic  force  in  the  direction  of  the  beam.  The 
acceleration  of  the  electrons  is  given  by  the  pendulum 


78 


equation,  and  thus  the  effect  of  tapering  the  undulator  is 
included  in  the  modified  pendulum  equation: 

C  =  8  +  |a|  cos(^  +  (f>)  (102. ) 

where  8  is  the  term  representing  the  new  'acceleration'  or 
'torque'  that  the  tapering  of  the  undulator  will  cause. 

To  determine  the  value  of  8  that  will  result  from 
reducing  the  undulator  wavelength  Xo,  start  with  equation 
(71.) 


V  =  i  =  L[(k  +  kj&  -  k]  (71.) 


and  calculate  the  effect  that  an  incremental  change  in  Ac 
will  have  on  the  phase  velocity  v  . 


v 


=  U 


k  + 


In 


lv 


\Pz 


8  = 
8  = 


.  /  2;r 

Av  =  L  (-  1)  — 

Ao  . 

A  „  M, 
Av  *  -2rtN — - 


PMo 


(103.) 
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where  the  last  line  included  L  =  nA0  and  fiz  ~  1 .  A  typical 


value  of  the  tapering  is  —  —  —0.2  .  If  N  —  100,  the  value 

for  the  acceleration  is  8  =  -200^r(-  0.2)  =  40?r  . 

The  value  of  8  achieved  by  a  change  in  the  magnetic 
field  strength  is  determined  similarly.  Start  with  equation 
(45.)  to  relate  the  undulator  parameter  Ku  to  fiz  by  using  a 
binomial  expansion: 


r'2( 1  +  K2)  =  1  -  Pi 


1- 

f  i  +  K) 

r2 

\  '  / 

1  +  K2 

2  y2 

(45.) 


(104.) 


Using  this  approximation  in  equation  (71.)  gives 


v  =  l  =  L[(k  +  k 0)PZ  -  k] 
v  =  L[k/?z  -  k  +  k0/?J 
V  =  L[k(/?z  —  l)  +  k0/?z] 


V  =  L 


1  - 


1  +  K2 


2  y 


\ 

ft  1  +  K2  ^ 

+  kn 

1 - 5-a- 

u 

) 

2r2 

^  '  /-I 

II 

-  k 

f  1  +  K2  'j 

+  k0 

l  J 

(105. ) 


80 


1  +  K2 

where  the  term  - —  «  1  was  dropped  from  the  second 

2  y 

parenthesis  term,  but  could  not  be  dropped  from  the  first 
parenthesis  term  because  it  was  the  only  remaining  piece. 

2n 

Continuing  by  making  the  substitutions  L  =  JSU0  ,  k  =  — ,  and 

A 

An 

using  equation  (29.)  for  the  value  of  —  the  value  for  the 

A 

acceleration  term  S  is 


=  Av  =  N/L 


2k  f  2k 


A  2  y 


A  A0  (  \  A 


1  +  K^2  y 

K2  \  AK 


1  L  o  AK,. 


i  +  K  K 


(106.) 


The  amount  of  tapering  that  can  be  applied  is  not 
unlimited.  The  general  rule  is 


|a|  >  |£|  <  4-yfa  >  2k 
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The  first  requirement  of  |a|  >  |£|  ensures  that  some  of 
the  electrons  will  be  trapped  in  the  ponderomotive  potential 
well,  known  as  particle  trapping.  If  a  is  smaller  than  5 
then  the  particles  will  all  be  accelerated  in  phase  space 
out  and  away  from  the  separatrix.  The  second  requirement 
ensures  that  the  tapering  is  strong  enough  to  overcome  the 
natural  untapered  deceleration,  and  the  last  requirement 
ensures  that  the  system  in  operating  in  a  strong  optical 
field  condition  where  tapering  is  needed.  [Ref.  29,  pg. 
157] 


G.  EFFICIENCY 

The  efficiency  T]  of  a  Free  Electron  Laser  is  defined  as 
the  fractional  loss  in  energy  of  the  electron  beam  during 
one  passage  through  the  undulator . 

_  (Ay)mc2  _  _(^r)  (107.) 

v  yt ic2  Yo 

where  (Ay)  is  the  average  energy  of  an  electron  at  the  end 
of  the  undulator. 

The  efficiency  of  the  FEL  can  be  increased  by  tapering 
the  undulator.  Assuming  X,  X0,  and  N  are  fixed,  the  effect 
of  tapering  ±  S  can  be  determined  starting  with  equation 
(105. ) 
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V  =  L 


(105.) 


f  1  +  Ku  1 

k  -  k 

o 

V  JJ 

Recalling  L  =  NX0  this  can  be  written  in  terms  of  the 
fixed  variables  X,X0,  and  N  as 


v 


tU0 


2  n  2 n 

fi  + 

X~T 

2  y2 

\  *  J- 

(108. ) 


Rearranging  terms,  the  equation  can  be  rewritten  as  a 
solution  of  y2 

v  = 

Y2  = 


For  a  tapered  undulator,  Ku  =  Kuo  +  AKur  .  Typically, 
AKU  «  Kuo  and  v  «  2m ,  so 


L  - 


M + 

2  y2X 

44  ±  Ku) 


1  - 


2m 


(109. ) 
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Y  = 


_  M  +  Kuo  +  2kuqAkut  +  (AKur)2) 


2M  1  - 


2m 


r  * 


.  A0(l  +  K20  +  2KuoAKur) 


1  + 


2/ZN 


2/^ 


1  - 


2;zN 


A0(l  +  K*0  +  2Ku0AKut]^1  + 


1  +  - 

v  2m  j 

-  •n 


2m 


(110.  ) 


2A 


where  higher  order  terms  involving  (AKur)2and  v2  have  been 

dropped.  Continuing  with  the  approximation,  equation  (110.) 
can  be  written  as 


y  2  _  Ap(l  +  K2j 

Y  2 A 


,  2k„.Ak  r 

V 

1  I  uo  u 

1  t  *7 - - — r 

1  +  - 

(!  +  Kuo)_ 

27ZN„ 

(111. ) 


Since  the  second  term  in  each  of  the  square  brackets  is 
assumed  small,  carrying  out  the  multiplication  yields 


*o(l  +  Kuo) 

"  .  2KuoAKur  V  ' 

1  +  -7 - - — r  *r 

2A 

(l  +  k  2J  2m] 

(112. ) 


At  t  =  0,  y  -  Y0  and  v  ~  vo> 


and  equation  (113.) 


becomes 
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(113.) 


2  =  ^ot1  +  Kuo) 

Yo  2/1 


1  +  ^ 


2m 


As  the  electrons  travel  down  the  undulator,  both  Yan<^ 
v  evolve  in  T,  and  away  from  their  initial  values  yo  and  vo  . 


Y  =  Yo  +  Ay 


(114.) 


and 


v  =  vQ  +  Av 


(115. 


Using  equation  (114.),  equation  (112.)  becomes 


y2  =  yI  +  2r0Ar  +  (Ar)2  * 


K  +  Kuo) 

22 


1  + 


2K  AK„r  vn  Av 

UP  u  _|_  O  . 

(i  +  k*J  2m  2m  j 


(116.) 


For  N  »  1,  (A y)2  «  yl ,  so  the  smaller  term  can  be 

dropped  from  equation  (116.).  Using  equation  (113.)  shows 
that  yl  can  be  cancelled  from  both  sides  leaving 


y2  =  ^y0Ay 


-  ^ot1  +  Kuo) 


2/1 


2KuoAKur  +  Av 

(l  +  kL)  2m 


(117. ) 
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Multiplying  both  sides  by 


yields 


Ay  /tc(l  +  K2J 

K  AK  T  Av 

UO  U  r 

7 — — - - — r  T 

Yo  2y2J 

_(l  +  K2J  4*nJ 

Now  using  equation  (106.)  to  define  the  taper  8 , 


8  = 


f 

-Am  ■ 


K, 


1  +  K,2 


\ 

/ 


AKu 

K 


equation  (118.)  becomes 


A  Y  A0(l  +  K20) 

8  +  Av" 

Yo  ~  2 ylA 

Am 

Using  the  resonance  condition  equation  (29.) 
y  =  yR ,  equation  (119.)  becomes 


Ay  _  yl  |~Av  -  <?" 
Yo  Yl  L  4^N  . 


For  N  »  1,  yo  ~  yR ,  and  equation  (120.)  simplifies  to 


(118. ) 


(106. ) 


(119. ) 


with 


(120. ) 
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(121. ) 


Ay  __  |~Av  -  8 

y0  L 


Equation  (107.)  can  be  written  as 


(Ay)  _  8  -  (Av) 
y  4tzN 


(122. ) 


where  Av  =  vf  -  vo ,  therefore 


„  =  l  Z  <v*>  + 

4^N 


(123.) 


where  vf  is  the  final  electron  energy,  and  vo  is  the  initial 
energy . 

Equation  (123.)  holds  for  both  positive,  8  >  0,  and 
negative,  <5  <  0 ,  tapering,  and  examples  of  the  change  in 
efficiency  for  various  taper  values  will  be  shown  in  the 
Chapter  V. 
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XV.  TJNAF  UPGRADE  SIMULATIONS 


In  order  to  move  from  damage  experiments  with  a  laser 
spot  size  of  a  few  millimeters  to  an  actual  weapon  with  a 
laser  spot  size  of  about  100  cm2,  the  average  power  of  the 
laser  must  be  increased.  The  current  average  power  at  the 
TJNAF  FEL  is  1.7  kW,  and  needs  to  be  increased  about  3 
orders  of  magnitude  to  several  megawatts  for  weapons 
applications . 

At  these  higher  MW  power  levels,  energy  efficiency 
becomes  important.  To  operate  more  efficiently,  the  TJNAF 
FEL  recirculates  the  electron  beam  after  giving  up  energy  to 
the  light  beam.  The  recirculated  electrons  are  timed  to 
arrive  back  at  the  entrance  to  the  RF  accelerator  180 
degrees  out  of  phase  with  the  alternating  electromagnetic 
field.  This  scheme  has  two  advantages.  First,  it  removes 
energy  from  the  electrons  so  that  they  can  be  sent  to  a 
final  beam  dump  without  risk  of  producing  neutron  radiation. 
This  requires  the  beam  energy  entering  the  dump  to  be  less 
than  10  MeV.  Second,  the  energy  removed  from  the  electrons 
can  be  used  to  boost  the  energy  of  the  RF  field.  A  diagram 
of  the  TJNAF  FEL  is  shown  in  Figure  10. 
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Recirculation 

adds 

at  least 

one 

new  constraint 

to 

the 

FEL  design.  Since 

the 

electrons 

are 

being 

returned 

to 

the 

entrance  to  the 

accelerator, 

they 

must 

undergo 

turns 

totaling  360  degrees.  The  turns  are  accomplished  through 
bending  magnets.  The  electron's  path  through  the  turns  is 
effected  by  its  momentum,  and  by  the  magnetic  force,  which 
is  proportional  to  the  magnetic  field  and  the  velocity  of 
the  electron.  The  magnets  are  designed  so  that  electrons 
moving  through  the  turn  at  the  correct  speed  will  travel 
through  the  center  of  the  beam  pipe  that  makes  up  the  turn. 
However,  electrons  with  higher  energy  will  be  traveling 
faster  than  this  design  speed,  and  above  a  certain  energy 
level  their  momentum  will  carry  them  into  the  outside  wall 
of  the  turn,  or  to  different  points  in  the  bending  magnet 
field,  or  to  different  points  at  the  entrance  of  the 
accelerator.  Similarly,  electrons  with  lower  energy  will 
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turn  too  quickly,  and  spread  the  beam.  These  effects  place 
limits  on  the  spread  in  energy  that  can  be  tolerated  for  a 
recirculated  beam,  and  will  have  to  be  considered  in  the 
upgrade  design. 

As  a  first  step  in  moving  from  1.7  kW  power  to  a 
megawatt  class  laser,  TJNAF  has  decided  to  upgrade  their 
current  design  to  produce  10  to  20  kW  average  power.  To 
support  this  goal,  simulations  were  conducted  to  determine 
gain,  power,  and  induced  electron  beam  energy  spread  that 
would  result  from  the  design  parameters  chosen  by  the  TJNAF 
personnel . 


A.  20  KW  SIMULATIONS 

The  first  sets  of  simulations  were  for  a  20  kW  design 
and  included  simulations  with  and  without  a  dispersive 
klystron  undulator  section.  As  shown  in  Figure  5  of  Chapter 
III,  an  undulator  can  have  a  dispersive  section  that  allows 
the  electrons  to  continue  transiting  the  optical  cavity 
without  any  interaction.  The  effect  of  this  dispersive 
section  is  to  allow  the  electrons  to  use  their  slightly 
altered  velocities  from  the  first  undulator  section  to  bunch 
electrons  prior  to  entering  a  second  undulator  section. 
This  results  in  a  significant  increase  in  gain  achieved  in 
the  second  undulator  section,  but  causes  a  reduction  in 
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saturated  power,  as  will  be  discussed  in  the  following 
sections . 

The  design  goals  are  to  achieve  20kW  of  average  power 
while  maintaining  a  final  energy  spread  of  6%  or  less  to 
ensure  that  recirculation  could  take  place.  In  order  to 
improve  gain  in  weak  optical  fields,  TJNAF's  20kW  design 
included  a  klystron  undulator.  The  simulations  varied  the 
dispersive  section  strength  D  of  this  klystron  undulator 
[Ref.  9],  and  varied  desynchronism  value  d  between  the 
optical  and  electron  pulse  in  order  to  compare  gain,  final 
power,  and  final  electron  energy  spread.  The  parameters 
used  in  the  simulations  are  shown  in  Table  4. 


Table  4.  20  kW  FEL  Parameters. 


PARAMETER 

Optical  Wavelength  A 

i 

*n 

o 

T— 1 

Undulator  Wavelength  A0 

20  cm 

Peak  Undulator  Magnetic  Field  B 

2kG 

Undulator  Periods  N 

24 

Undulator  Length  L 

6  m 

Undulator  Parameter  K 

0.79  1 

Dispersive  Strength  D 

__  0,  1,  2  or  3 

Electron  Energy  ymc2 

200  MeV 

Electron  Beam  Radius  rb 

0.5  mm 

Electron  Pulse  Length  le 

2  ps 

Initial  Electron  Phase  Velocity  Spread 

Standard  Deviation  crG 

0.3 

Electron  Injection  Angle  Standard  Deviation  as 

0.15 

Peak  Current  I 

200  A 

Average  Current 

8  mA 
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2 

Cavity  losses  (Q  =  10) 

10% 

Peak  Optical  Power  within  cavity 

40  GW 

Repetition  Rate 

74.85  MHz 

Duty  Cycle 

4x1 0'5 

1 .  Gain 

Simulations  were  conducted  for  values  of  dispersive 
strength  D  =  0,  1,  2,  and  3,  with  D=0  representing  the 

conventional  undulator  with  no  dispersive  section.  The 
desynchronism  between  electron  pulse  and  optical  pulse  was 
also  varied  from  d  =  0.0  to  0.30  in  steps  of  0.02. 

The  output  from  a  typical  simulation  is  shown  in  Figure 
11.  The  variables  used  in  this  particular  simulation  are 
shown  across  the  top  of  the  figure  and  are:  dimensionless 

current  j  =2,  the  electron  pulse  length  az  -  =  ^ ' 


angular  electron  velocity  spread  <Jg 


0.15 


where  A6  is  the  electron  beam  angular  spread,  electron 
phase  velocity  energy  spread  <rG  =  4^N  ^ =  0.3,  where 
is  the  fractional  energy  spread,  desynchronism 


d  =  =  0.03 ,  resonator  cavity  quality  factor  of  Q  =  106 

to  represent  a  cavity  with  practically  no  losses,  undulator 
periods  N  =  24,  and  dispersive  section  strength  D  =  0.  The 
high  Q  factor  is  chosen  to  determine  the  gain  produced  from 


93 


the  beam  interaction  alone,  ignoring  any  loses  in  the 
resonant  cavity. 


-2  z  2  -13  V  13  -13  V  13 


-2  z  2  -13  V  13  0  n  40 


Figure  11.  Simulation  Results  Showing  Gain. 


Figure  11  shows  9  sections  arranged  in  a  3  x  3  grid 
pattern.  The  lower  left  corner  indicates  the  position  of 
the  parabolic  electron  pulse  shape.  The  electron  pulse 
slips  back  by  the  slippage  distance  NX  relative  to  the  light 
pulse  as  both  pulses  move  through  the  undulator  from  x=  0 
to  T  =1 .  The  middle-left  plot  shows  the  evolution  of  the 
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optical  field  |a(z,  n)|  as  it  grows  over  n  =  40  passes  through 
the  undulator,  with  the  final  optical  pulse  shape  shown  in 
the  top- left  window. 

The  center  graph  shows  the  evolution  of  the  optical 
power  spectrum,  P(v)  over  n  =  40  passes  through  the 
undulator.  The  top-center  plot  shows  the  final  optical 
power  spectrum  after  n  =  40  passes.  The  bottom-middle  plot 
is  the  weak-f ield  gain  spectrum,  G(v),  showing  a  maximum 

gain  occurring  at  v  =  2.6. 

The  lower-right  window  shows  the  gain  of  the  optical 
pulse  G  (n)  .  The  middle-right  plot  shows  the  evolution  of 
the  electron  spectrum  f  (v)  over  n  =  40  passes,  while  the 
top-right  plot  shows  the  final  spectrum  after  n  =  40  passes 
through  the  undulator. 

The  gain  results  from  many  simulations  at  different 
values  of  desynchronism  d  is  plotted  in  Figure  12.  The 
conventional  undulator  with  D  =  0  achieved  a  maximum  gain  of 
G  =  0.2  at  a  desynchronism  value  d  =  0.03.  The  maximum  gain 
of  G  =  0.67  was  achieved  for  the  klystron  undulator  with 
dispersive  strength  D  =  2 .  This  was  3  times  higher  than  the 
gain  of  the  conventional  undulator  with  D  =  0.  However,  the 
higher  gain  achieved  using  a  klystron  undulator  comes  at  the 
expense  of  a  lower  power  at  saturation,  as  described  in  the 
next  section. 
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0.8 


Figure  12.  Steady  State  Gains 


2 .  Power 

A  typical  output  showing  the  power  evolution  of  the 
optical  beam  is  shown  in  Figure  13 .  All  the  variables  used 
in  this  simulation  were  the  same  as  used  in  the  gain 
simulation  of  Figure  11,  except  that  the  simulation  has  been 
allowed  to  run  to  n  =  1200  passes,  and  the  quality  factor 
has  been  reduced  to  Q  =  10.  This  Q  represents  the  actual 
TJNAF  FEL  cavity  loss  of  1/Q  =  10%.  In  Figure  13,  the 
optical  power  as  a  function  of  n  passes  is  plotted  in  the 
lower-right  corner,  while  all  the  other  windows  of  the  plot 
remain  the  same  as  described  in  Figure  11. 
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-5  2  S  -47  V  47  0  n  1200 

Figure  13.  Simulation  Results  for  Conventional  Undulator. 

Figure  14  shows  the  results  of  many  simulations 
plotting  the  final  average  power  as  a  function  of  dispersive 
strength  D  and  desynchronism  value  d.  Maximum  power  was 
achieved  with  the  conventional  undulator  and  dispersive 
strength  D  =  0.  The  power  level  achieved  was  14  times 

higher  than  the  closest  klystron  undulator  which  used 
dispersive  strength  D  =  1.  Klystrons  with  dispersive 

strengths  D  =  2  and  D  =  3  resulted  in  even  lower  power 
levels . 

The  results  shown  in  Figure  14  seem  to  indicate  that 
the  conventional  undulator  is  the  obvious  design  to  choose 
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in  order  to  achieve  maximum  power.  However,  strong  fields 
in  the  conventional  undulator  cause  the  trapped-particle 
instability.  This  causes  the  electron  beam  energy  spread  to 
move  outside  the  6%  design  tolerance. 


Figure  14.  Final  Power  vs.  Desynchronism 

3 .  Energy  Spread 

Figure  15  shows  the  final  induced  electron  beam  energy 
spread  produced  by  the  4  different  undulators  for 
desynchronism  values  d  =  0.0  to  0.15.  As  can  be  seen  in 
this  figure,  the  conventional  undulator  with  D  =  0  has  a 
peak  energy  spread  of  8.6%  occurring  at  d  =  0.007.  The 
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energy  spread  decreases  with  increasing  de synchronism,  and 
falls  below  the  desired  6%  level  at  a  desynchronism  value  of 
d  =  0.06,  and  remains  below  6%  for  all  larger  values  of  d. 

For  the  klystron  undulator  with  dispersive  strength 
D  =  1,  the  maximum  energy  spread  is  5%,  which  occurs  at 
desynchronism  values  from  d  =  0 . 01  to  0 . 06 .  It  also  shows  a 
decrease  in  the  energy  spread  with  an  increase  in 
desynchronism  above  d  =  0.06.  The  klystrons  with  D  =  2  and 
D  =  3  have  an  energy  spread  of  about  2%  for  all  values  of 
de  synchr on i sm . 


Fractional  Spread  in  Energy  vs.  Desynchronism. 


Figure  15.  Energy  Spread  vs.  Desynchronism 
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4.  Conclusions  of  20kW  Simulations 

As  expected,  the  maximum  gain  was  achieved  with  a 
klystron  undulator.  Using  a  klystron  undulator  with 
dispersive  strength  D  =  2  produced  a  steady  state  gain  of 
67%.  However,  this  increased  gain  leads  to  a  tremendous 
reduction  in  final  power.  The  D  =  2  klystron  design's  final 
power  was  approximately  15  times  smaller  than  that  of  the 
conventional  undulator. 

The  conventional  undulator  achieved  the  highest  final 
power,  which  occurred  at  a  desynchronism  of  d  =  0.007. 
However,  at  this  desynchronism,  the  induced  energy  spread 
was  8.6%,  which  is  outside  the  limits  imposed  by 
recirculation.  Increasing  the  desynchronism  to  d  =  0.03 
reduces  the  energy  spread  down  to  the  desired  6%,  while 
achieving  a  final  output  power  that  is  still  9  times  higher 
than  the  best  klystron  undulator  design.  Simulations  using 
a  conventional  undulator  with  a  desynchronism  value  of  d  = 
0.03  produced  a  final  power  of  19.2  kW,  which  is  close  to 
the  20kW  objective,  and  corresponds  to  a  final  gain  of  19%. 

Results  for  the  20  kW  simulations  were  presented  at  the 
20th  International  Free  Electron  Conference  in  Williamsburg, 
VA,  and  were  published  in  Nuclear  Instruments  and 
Measurements.  [Ref.  9] 
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B.  10  KW  SIMULATIONS 

Recently,  simulations  for  a  10  kW  upgrade  design  were 
conducted.  The  design  parameters  are  listed  in  Table  5. 


Table  5.  10  kW  FEL  Parameters 


PARAMETER 

Optical  Wavelength  A 

5.8  pm 

Undulator  Wavelength  A0 

20  cm 

Undulator  Periods  N 

24 

Undulator  Length  L 

4.8  m 

Undulator  Parameter  K 

1, 2,  3,  or  4 

Dispersive  Strength  D 

0, 0.25, 0.5, 0.75  or  1 

Electron  Energy  ymc2 

150  MeV 

Electron  Beam  Radius  rb 

0.35  mm 

Electron  Pulse  Length  az 

0.5  ps 

Initial  Electron  Phase  Velocity  Spread 

Standard  Deviation  aG 

0.5 

Peak  Current  1 

270  A 

Average  Current 

10  mA 

Dimensionless  Current  Density  j 

2.9  -  3.7 

Cavity  losses 

15% 

Peak  Optical  Power  within  cavity 

1.5  MW 

Repetition  Rate 

74.85  MHz 

Duty  Cycle 

3.75x10  s  1 

In  addition  to  conducting  simulations  for  various 
values  of  dispersive  strength  and  desynchronism,  these 
simulations  also  varied  the  value  of  the  undulator  parameter 
K  from  1  to  4  in  steps  of  unity.  Table  6  shows  the 
parameters  affected  when  the  undulator  parameter  is  varied: 
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X =  the  optical  wavelength,  =  the  dimensionless  pulse 
length,  j  is  the  dimensionless  current  density,  and  F  is  the 
filling  factor  which  is  the  ratio  of  the  area  of  the 
electron  beam  to  the  area  of  the  optical  beam  averaged  over 
the  interaction  region. 


Table  6.  Dimensionless  Variables  vs.  K. 


K=1 

K=2 

K=3 

K=4 

2  f3  JXm 

5.8 

12  pm 

20  pm 

2.7 

1.1 

0.54 

0.32 

jF 

2.9 

3.7 

3.8 

3.8 

F 

0.06 

0.02 

0.01 

7xl0"3 

1 .  Power 

Simulations  were  carried  our  for  desynchronism  values  d 
ranging  from  0  to  0.3.  The  dispersive  strengths  for  these 
simulations  were  varied  from  D  =  0.0  to  D  =  1.0  in  steps  of 
0.25.  Figure  16  shows  the  output  for  the  conventional 
undulator  design  D  =  0,  desynchronism  d  =  0.025,  and 
undulator  parameter  K  =  2 . 
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Figure  16.  Power  Simulation  for  Klystron  Undulator 


Figure  17  shows  the  maximum  power  achieved  for  various 
values  of  K  and  D  =  0  (no  klystron)  .  The  maximum  power 
was  achieved  for  the  K  =  2  undulator  parameter,  with  a  final 
power  of  14.5  kW.  This  was  also  the  maximum  power  achieved 
out  of  all  360  simulations  conducted  for  the  various 
undulator  parameters  K,  dispersive  strengths  D,  and 
desynchronism  values  d.  Figure  17  also  shows  that  the  K  =  1 
simulations  achieved  similar  final  power  levels  at  varying 
desynchronism  values.  The  energy  spread  for  these 
simulations  was  approximately  5.5%,  which  was  within  the  6% 
specification  for  recirculation.  No  power  is  achieved  for 
the  K  =  3  and  4  simulations  due  to  the  smaller  filling 
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factor  F  and  the  shorter  dimensionless  electron  pulse  length 


oz.  Note  that  the  actual  electron  pulse  length  le 
fixed,  but  the  slippage  distance  NX  is  larger. 


Figure  17.  Final  Power,  Conventional  Undulator 


remains 
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15 


1 


D=0.25 


Figure  18.  Final  Power,  Weak  Klystron  Undulator 


Figure  18  shows  the  final  power  after  introduction  of  a 
weak  klystron  of  D  =  0.25.  The  figure  shows  that  the 
klystron  has  added  enough  gain  to  produce  power  for  the  K  = 
3  simulations,  but  at  a  final  power  level  of  only  7.5  kW. 
The  maximum  power  achieved  using  the  K  =  1  and  2  values  has 
dropped  from  14.5  kW  to  approximately  10  kW.  The  energy 
spread  for  these  simulations  was  reduced  to  about  4%  because 
of  the  lower  power. 

With  stronger  klystron  values  of  D  =  0.5  and  up  to  D  = 
1,  all  four  undulator  parameter  values  produced  some  output 
power,  but  less  than  the  desired  10  kW.  The  maximum  average 
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power  achieved  for  the  higher  D  values  was  for  K  =  3,  D  = 
0.50,  which  produced  a  peak  power  of  6.7  kW. 

2.  10  kW  Simulation  Conclusions 

From  the  360  simulations  conducted  for  the  10  kW 
design,  it  appears  that  the  conventional  undulator  D  =  0, 
and  the  weak  klystron  undulator  with  D  =  0.25,  will  meet 
the  desired  goals.  However,  both  undulator  designs  have 
marginal  areas  of  operation.  For  the  weak  klystron 
undulator  of  D  =  0.25,  the  energy  spread  was  very  good  at 
4%,  but  the  final  output  power  was  only  slightly  higher 
than  the  desired  10  kW.  For  the  D  =  0  undulator  (no 

klystron),  the  output  power  was  much  higher  at  14.5  kW,  but 
the  energy  spread  was  about  5.5%,  which  is  close  to  the  6% 
limit . 

Results  for  the  10  kW  simulations  were  presented  at 
the  22nd  International  Free  Electron  Laser  Conference  in 
Raleigh,  NC,  and  will  be  published  in  Nuclear  Instruments 
and  Measurements.  [Ref.  7] 
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V. 


TAPERING 


It  is  only  near  resonance  that  the  coupling  between  the 
laser  and  electron  beams  is  significant.  As  the  electrons 
lose  energy  to  the  light  beam,  they  move  away  from 
resonance.  In  order  to  maintain  or  re-establish  the 
resonance  condition,  the  undulator  properties  can  be  tapered 
so  that  as  the  electrons  lose  energy,  the  resonant  energy 
also  decreases.  The  theory  of  tapering  by  either  changing 
the  undulator  wavelength  Xa,  or  changing  the  undulator 
magnetic  field  strength  is  discussed  in  section  F  of  Chapter 
III. 

Without  tapering,  the  electrons  begin  to  bunch  early  in 
the  growing  optical  field,  and  eventually  move  to  points  in 
phase  space  where  they  are  no  longer  losing  energy  but  are 
gaining  energy  from  the  light  beam.  When  the  optical  field 
reaches  this  level  the  FEL  is  said  to  have  reached 
saturation.  The  tapered  -undulator  increases  gain  and 
efficiency  in  strong  optical  fields  and  extends  the  usual 
saturation  limit  to  stronger  optical  fields  by  maintaining 
the  resonance  condition.  As  electrons  lose  energy  ^mc2'  the 


resonance  condition. 


would  indicate  that 
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either  reducing  the  undulator  wavelength,  X0,  or  the 

undulator  parameter,  K,  would  restore  resonance. 

Since  the  purpose  of  tapering  is  to  reduce  the 
resonance  energy  to  match  the  reduced  energy  of  the 
electrons,  it  would  seem  that  there  would  be  no  benefit  in 
using  an  inverse  taper  where  the  undulator  wavelength, 
is  increased,  or  the  undulator  magnetic  field  strength,  K, 
is  increased.  However,  simulations  show  that  there  are 
several  possible  advantages  to  using  the  inverse  taper. 

The  electron  dynamics  in  a  tapered  undulator  were 
described  by  equation  (102.)  of  chapter  III  as 

oo 

£  =  8  +  |a|  cos(^  +  </> )  (102.) 

where  f  is  the  electrons  phase  in  the  ponderomotive  wave,  (/> 
is  the  optical  phase,  a  =  |a|ei0  is  the  complex  dimensionless 

optical  field  strength,  and  8  is  the  acceleration  torque 
term  resulting  from  tapering.  Simulations  using  both 

positive  and  negative  values  of  8,  corresponding  to  positive 
and  negative  taper  of  K  or  Xo,  were  carried  out  and  the 
results  will  be  discussed  in  following  sections. 
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A.  EFFICIENCY  AND  ENERGY  SPREAD 

When  designing  a  high  power  recirculating  free  electron 
laser,  both  efficiency  and  induced  energy  spread  are 
important.  As  defined  in  equation  (107.)  of  Chapter  III, 

efficiency  7]  is 


(107. ) 


and  is  typically  highest  in  strong  fields.  However,  strong 
fields  have  the  adverse  effect  of  inducing  a  larger  electron 
energy  spread  along  the  undulator.  For  our  purposes,  energy 
spread  is  defined  as  the  full  spread  in  phase  velocity, 


The  difficulty  for  a  recirculating  FEL  operating  in 
strong  fields  arises  because  some  electrons  gain  energy 
while  some  lose  energy.  This  makes  the  spread  in  energy 
large  but  does  not  necessarily  make  the  efficiency  large. 

A  second  difficulty  arises  when  considering  high  power 
FELs.  The  small  percentage  of  electrons  that  may  lie  in  the 
outer  tails  of  the  energy  distribution  cannot  be  ignored, 
because  for  high  power  FELs  these  tails  represent  a  lot  of 
current.  When  recirculating  high  current  (for  high  power). 
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this  small  fraction  of  the  electron  beam  in  distribution  can 
cause  serious  damage  to  the  vacuum  beam  pipe,  as  well  as 
radiation  problems  as  the  electrons  collide  with  the  walls 
of  the  recirculating  path.  Consider  just  0.01%  of  a  1A 
(average  current)  beam  at  100  MeV  energy  focussed  to  a  1  mm2 
spot;  the  irradiance  is  1  MW/cm2.  The  damage  would  occur 
rapidly.  Therefore,  when  discussing  energy  spread  for 
high  current  recirculating  devices,  the  electrons  in  the 
tails  of  the  distribution  must  be  included  as  in  equation 
(124. )  . 

B.  NO  TAPER 

In  the  conventional  undulator  with  weak  initial  optical 
field  ae  <  n,  the  gain  spectrum  is  shown  in  Figure  19.  It 
is  symmetric  about  resonance  with  a  peak  gain  of  27%  for  j  = 
2,  occurring  at  an  initial  phase  velocity  vD  =  2.6.  The 
gain  spectrum  width  is  vs  ~  n ,  and  the  gain  is  zero  at 

resonance  vo=  0 . 
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Peak  Gain  of 
27%  at 
v0  =  2.6 


As  the  optical  field  grows,  the  gain  spectrum  changes 
in  shape.  Figure  20  shows  the  spectrum  as  the  optical  field 
grows  from  ao  =  small  to  saturation  at  ao  =  40. 


ln(l-K5) 


Figure  20.  Gain  vs.  v0  in  Increasing  Optical  Field  a0 


The  spectrum  is  still  symmetric  about  resonance,  but 
the  peak  gain  has  been  reduced  to  approximately  1%  at  an 
initial  phase  velocity  of  VQ  =  10.  At  saturation,  the  FEL 
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oscillator  wavelength  will  shift  to  follow  the  point  of 
maximum  gain,  and  in  this  case  will  shift  from  vo  =  2.6  to 

v0  =10,  which  corresponds  to  a  shift  of  about  5%  in  optical 
wavelength  for  N  =  25. 

In  the  electron  pendulum  phase- space  determined  by 
Equation  (102.),  saturation  in  strong  fields  occurs  when 
electrons  become  trapped  in  closed  orbits  as  shown  in  Figure 
21. 


— Jt/2  £  3n/2  0  T  1 

Figure  21.  Phase  Space  for  Conventional  Undulator. 
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When  the  pendulum  separatrix  height,  4-\/a  ,  becomes 
greater  than  the  phase  velocity  for  peak  gain  in  weak 
fields,  vo  =  2.6,  an  increasingly  larger  fraction  of  the 

electron  beam  becomes  trapped  in  the  closed  orbit  region  of 
phase  space.  At  saturation,  the  FEL  interaction  induces  a 
spread  in  phase  velocities  which  is  slightly  larger  than  the 

height  of  the  separatrix,  Vs  ~  4-v/a  *8 71,  or  about  10%  for  N 
=  25.  This  spread  is  sufficient  to  effect  the  recirculation 
of  the  electron  beam. 

C.  POSITIVE  TAPERING 

When  someone  speaks  of  tapering  the  undulator,  they 
usually  mean  positive  tapering:  either  the  undulator  period 
X0  is  decreased  along  the  undulator,  or  the  magnetic  field 
strength  is  decreased  along  the  undulator,  or  both.  With  a 
modest  taper,  the  maximum  weak  field  gain  is  shifted  and 

occurs  at  initial  phase  velocity  vG  «  2.6  —  ^  •  As  can  ke 

seen  by  comparing  Figure  20  with  Figure  22,  the  weak-field 
gain  is  reduced  for  a  tapered  undulator.  However,  strong 
field  gain  is  increased.  For  a  taper  of  6n,  the  weak  field 
maximum  gain  is  shifted  to  vo  ~  —7  .  In  strong  fields  of 
ao  =  40,  the  gain  spectrum  has  become  flattened  and 
distorted,  and  the  maximum  gain  has  shifted  to  vo  ~  0 .  In 
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strong  fields  the  gain  has  increased  from  G  =  1%  in  the 


untapered  undulator,  to  G  «  2%  with  taper. 


Peak  Weak  Field  Gain  Shifts  left  by  5/2 
and  is  reduced  from  the  Untapered  Gain. 


Figure  22.  Gain  vs.  v„  in  Increasing  Optical  Field  a0 
for  a  Tapered  Undulator. 


The  separatrix  surrounding  the  closed-orbit  paths  in 


phase  space  is  given  by  [Ref.  29,  pg.  157] 


VS  =  2S(Cs  -  Co)  +  2a(sin (£.)  -  sin(C0))  (125.) 


where 


C0  =  In  -  cos  '(-  5  /|a|) 


(126. ) 
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In  the  tapered  undulator,  the  phase  paths  are  modified 
by  the  acceleration  term  8  in  the  pendulum  equation.  For 
8  >  a,  the  cos'1  expression  has  no  solution,  and  there  is  no 

separatrix  surrounding  closed  orbit  paths. 

The  change  in  the  separatrix  can  be  calculated  starting 
with  equation  (125.).  With  a  tapered  undulator,  the 
separatrix  height  is  reduced,  and  has  been  shown  by  Massey 
[Ref.  30]  to  be 


h  =  4*\/a^|q|cos”1  (-  |q|)+  -y/l  -  q2  j 


(127. ) 


where  q  =  ^  . 


When  the  taper  is  too  large,  such  that 


8>a0,  the  electrons  phase  space  has  no  closed  orbits. 
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8  >a0 

Accelerating  electrons  never 
enter  closed  orbits. 


Time  evolution  is  shown 
as  electrons  change  from 
light 
to  dark. 


Figure  23.  Phase  Space  Showing  8  >  a0. 


In  sufficiently  strong  fields,  ao  >  S,  some  of  the 

electrons  can  be  trapped  in  the  closed  orbits  of  the 
pendulum  phase  space.  This  is  illustrated  in  Figure  24.  In 
stronger  fields,  ao  »  S,  a  large  fraction  of  the  electrons 
can  be  trapped  near  resonance  around  £  ~  n  ,  and  continue  to 
lose  energy  to  the  optical  field,  as  is  illustrated  in 
Figure  25. 
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5»a0 
Most  electrons 
enter  closed  orbits. 


Time  evolution 
shown 

as  electrons 
change  from  lig)tit 
to  dark. 

“20 


Gain 


0.1 


0 

0.05 


J.l  j.l  UU  U.LU-U-U4. 

3n/2  Ox  1 


Figure  25.  Phase  Space  Showing  a0  »  5. 


As  shown  in  Figure  26  below,  electrons  starting  near 
the  phase  £  ~  0  are  accelerated  by  both  the  torque  S  and  the 

strong  optical  field  a,  while  the  electrons  starting  near 
£  =  n  see  the  acceleration  term  8  and  a  roughly  cancel, 

leaving  them  trapped  in  closed  orbits.  Tapering  is 
effective  because  the  electrons  near  the  maximum  gain  point 
of  £  ~  n  are  trapped,  while  the  electrons  near  the  phase 

point  for  absorption  £  ~  0 ,  are  taken  away  from  resonance 
and  eventually  stop  interacting.  It  is  this  imbalance 
between  the  gain  and  loss  mechanisms  that  leads  to  gain. 
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Torque  term  5  and 
optical  Held  aG  cancel 
together  for  electrons 

Torque  term  5  and  starting  near  £0  =  n- 

optical  field  a0add 
together  for  electrons 
starting  near  £o=0. 


Figure  26.  Separatrix  for  Positive  Tapered  Undulator. 


As  can  be  seen  in  Figure  27,  which  shows  only  the  final 
position  of  each  electron,  the  electrons  appear  to  be  have 
an  energy  spread  with  about  half  of  the  electrons  trapped 
near  resonance  v  ~  0  ,  and  the  other  half  untrapped  at 

V  ~  8  +  2-s/a  .  In  a  strong  field  of  a  =  40  and  with  a  taper 
rate  of  S  =  5  n,  the  induced  spread  approximately 

VS  =  <?  +  2V^  15  +  2740=30.  Due  to  this  large  energy  spread, 

the  tapered  undulator  does  not  appear  to  be  desirable  for 
recirculating  the  electron  beam. 
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-n/2  C,  3n/2  0  t  1 


Figure  27.  Phase  Space  for  Positive  Tapered  Undulator. 


D.  NEGATIVE  TAPER 

In  the  negative  or  inverse  tapered  undulator,  the 
acceleration  term  8  <  0.  Initial  investigations  and 

simulations  into  the  use  of  negative  taper  indicated  some 
favorable  results  could  be  obtained  in  the  areas  of  strong 
field  gain  and  energy  spread. 

Figure  28  shows  the  gain  spectrum  for  an  FEL  with 
inverse  taper  rate  8  -  -6n  over  N  =  24  periods  with  a 
current  density  of  j  =  2 .  These  parameters  are  descriptive 
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of  the  TJNAF  FEL  with  recirculation  of  the  electron  beam. 
In  weak  fields,  a  <  n ,  the  gain  spectrum  peak  is  shifted 

above  resonance  to  vo  ~  2.6  —  ^  =  12,  with  a  peak  gain  of 

G  =  12%.  This  gain  is  comparable  to  the  positive  taper  case 
shown  previously  in  Figure  22  above.  Both  positive  and 
negative  taper  have  the  effect  of  shifting  the  gain  spectrum 
by  vo  ~  2 .6  -  ^ ,  and  reducing  the  peak  gain  available  by 

about  50%  when  8  =  ±6n  . 

In  stronger  fields,  the  gain  spectrum  distorts  and 
flattens,  but  in  a  manner  that  differs  from  both  the  non- 
tapered  and  positive  tapered  undulators.  As  shown  in  Figure 
28,  the  peak  in  the  gain  spectrum  tends  to  stay  near  the 
same  value  vo  with  only  a  slight  shift  from  vo  =  12  to 

vo  =  13.5  as  the  field  grows  from  weak  to  strong  values. 

This  means  that  in  a  FEL  with  inverse  taper,  the  expected 
wavelength  shifts  only  slightly  as  the  laser  approaches 
saturation.  The  peak  gain  in  strong  fields  does  decrease, 
but  only  down  to  3%,  which  provides  a  significant  advantage 
over  the  untapered  undulator,  and  even  betters  the  positive 
tapered  undulator.  This  result  was  unexpected  since  the 
tapered  undulator  has  long  been  thought  to  provide  the  best 
extension  of  the  FEL  to  strong  fields. 
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Strong  Field  Gain  tends  to  remain  near 
weak  Field  gain  maximum  v0.  Strong  Field 
In  (1+G)  decreases  but  to  a  level  significantly 
.  higher  than  the  conventional  and  positive 

]  tapered  undulators. 


Figure  28.  Gain  vs.  v0  in  Increasing  Optical  Field  a0 
for  a  Inversely  Tapered  Undulator. 


The  electron  evolution  in  phase  space  with  inverse 
taper  is  also  quite  different  than  the  other  two  undulator 
types.  The  separatrix  shown  in  Figure  29  is  still  given  by 
equation  (124.)  above,  but  with  8  <  0.  Figure  29  is  shown 

for  a  taper  rate  of  8  =  -6?r  ,  N  =  24  periods,  and  a  current 
density  of  j  =2. 
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Electrons  start  at  v0  =  15. 


Electrons  cannot 
enter  closed  seperatrix, 

and  must  bunch  as _ 

they  move  around 
separatrix  to  low 
energy  level,  resulting 
in  coherent  radiation. 

-tt/2  C  331/2  0  T  1 

Figure  29.  Phase  Space  for  an  Inverse  Tapered  Undulator. 

The  electrons  started  at  vo  =  15  ,  which  provides  the 

peak  gain  in  strong  fields  for  this  taper.  All  of  the 
electrons  start  in  open  phase  space  paths  far  from  resonance 
and  cannot  enter  the  closed-orbit  region  inside  the 
separatrix.  Following  the  open-orbit  paths,  it  is  clear 
that  the  electrons  must  bunch  as  they  lose  energy  to  travel 
around  the  closed-orbit  region.  The  bunch  continues  to 
follow  open  orbit  paths  down  to  a  final  phase-space  site 

around  V  ~  -20  where  they  have  lost  considerable  energy. 

More  importantly,  the  induced  average  energy  spread  in  the 
final  electron  distribution  appears  to  be  relatively  small. 
Most  of  the  final  electrons  are  located  in  a  spread  of 
vs  ~  15  around  v  =  -16,  with  a  tail  extending  back  to 
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resonance.  This  corresponds  to  only  a  4%  average  energy 
spread.  This  is  significantly  reduced  compared  to  the 
tapered  undulator  spread  of  13%  and  the  periodic  undulator 
spread  of  10%. 

These  results  are  surprising  and  appear  promising  for  a 
recirculating  FEL.  For  high  power  recirculating  FELs,  it 
is  the  full  energy  spread  that  is  important.  The  electrons 
at  the  extreme  edges  of  the  distribution  must  be  considered. 
Figure  30  provides  a  phase  space  view  of  a  undulator  with  a 
taper  of  8  =  -5n  ,  showing  only  the  final  position  of  the 
electrons.  Most  of  the  electrons  are  bunched  around  V  =  - 
20,  there  are  still  some  electrons  at  the  top  of  the 
separatrix  around  v  =  8 .  When  considering  the  full  spread 
of  electrons,  a  comparison  of  Figure  21  and  Figure  27  with 
that  of  Figure  30  shows  that  the  inverse  tapered  undulator' s 
full  spread  is  vs  ~  28 .  This  is  only  a  modest  improvement 
over  the  untapered  and  positive  tapered  undulator  with  a 
full  spread  of  vs  =  30  . 
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Figure  30.  Phase  Space  for  a  Negative  Tapered  Undulator. 


E.  CONCLUSIONS  OF  THE  CONVENTIONAL  AND  TAPERED  UNDULATORS 

The  inverse  tapered  undulator  appears  to  offer 
advantages  over  both  the  periodic  and  tapered  undulators . 
The  advantages  include  larger  gain  in  strong  fields,  smaller 
wavelength  shift  at  saturation,  and  a  smaller  induced 
average  energy  spread. 

The  increase  in  gain  over  the  untapered  undulator  in 
strong  fields  appears  to  be  about  3  to  1,  while  the  increase 
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in  gain  over  the  negative  tapered  undulator  is  about  1.5  to 

1. 

While  both  the  untapered  and  positive  tapered 
undulators  have  large  shifts  in  resonant  wavelength  as  the 
optical  field  grows  to  saturation,  the  inversely  tapered 
undulator  experiences  only  a  minor  shift.  The  result  for 
inverse  taper  is  only  a  slight  shift  in  the  lasing 
wavelength  as  the  optical  field  grows  from  start-up  to 
saturation. 

The  induced  full  width  half  maximum  energy  spread  is 
greatly  reduced  for  the  inverse  tapered  undulator.  However, 
the  full  width  energy  spread  is  only  slightly  improved  by 
using  the  inverse  taper  design. 

Initial  results  from  the  inverse  taper  simulations  were 
presented  at  the  21st  International  Free  Electron  Laser 
conference,  and  published  in  Nuclear  Instruments  and 
Measurements.  [Ref.  8]  Additional  results  were  presented  at 
the  22nd  International  Free  Electron  Laser  conference,  and 
will  be  published  in  Ref.  7. 


126 


LIST  OF  SYMBOLS 


a  Dimensionless  optical  field 

o 

a  Dimensionless  derivative  of  optical  field 
ao  Dimensionless  optical  field  at  start  of  undulator 
A  Area 

Bx/y/  z  Magnetic  Field  in  x,  y,  or  z  direction 
Bu  Undulator  Magnetic  Field 

P  Fractional  speed  (velocity/speed  of  light) 
c  Speed  of  light  in  vacuum 

C  Specific  Heat 

d  desynchronism  value 

D  Dispersive  strength 

Dl  Thermal  Diffusion  length 

8  Dimensionless  taper 

83  Three  dimensional  Dirac  delta  function 

e  Electron  charge 

E  Energy 

Es  Electric  Field 

8  Unit  vector  symbol 

F  Filling  factor 

O  Irradiance  or  Flux  density 

<I>o  Required  Flux  Density 

0  Initial  phase  of  optical  wave 
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G  Optical  field  gain 

Y  Lorentz  factor 

Y0  Initial  Lorentz  factor 
Yr  Resonance  condition  Lorentz  factor 
AHm  Latent  heat  of  melting 
AHv  Latent  heat  of  vaporization 
r]  Efficiency 

j  Dimensionless  current 
J  Current 

k  Optical  wave  vector 

ko  Undulator  wave  vector 
K  Thermal  conductivity 
Ku  Undulator  parameter 

K  Thermal  diffusivity 

L  Length 

X  Optical  wavelength 

Xc  Undulator  wavelength 
m  Mass 

n  Number  of  passes  through  the  undulator 
N  Undulator  Periods 

v  Phase  velocity 

vf  Final  Phase  velocity 
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vo  Initial  Phase  velocity 

vs  Energy  Spread 

P  Power 

Q  Resonator  quality  factor 

r2  Distance  in  the  'i'  direction 

R  Radius 

RiH  Radius  of  the  hole  for  the  'ith'  run 

RiT  Radius  of  the  entire  damaged  area  for  the  'ith' 

p  Mass  density 

pe  Electron  density 

S  Resonator  cavity  length 

oG  Electron  phase  velocity  spread 

Ce  Electron  angular  spread 

oz  Electron  pulse  length 

T  Temperature 

Tm  Melting  temperature 

Tc  Initial  temperature 

Tv  Vaporization  temperature 

Tp  Pulse  length 

T  Dimensionless  variable  for  length 
v  velocity 

V  Volume 

ViH  Volume  of  the  hole  for  the  'ith'  run 

ViT  Volume  of  the  entire  damaged  area  for  the  'ith' 


run 


run 
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(0 


Radian  frequency 


VP  Phase  of  optical  wave 

y  Thickness  of  material 

z  Depth  of  hole 

£  Electron  Phase 
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